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We  undertake  a  systematic  study  of  the  adiabatic  dynamics  of 

nonlinear  excitations  within  the  Su,  Schrieffer  and  Heeger  (SSH)  model 

for  polyacetylene.   The  SSH  model  is  a  tight-binding  electron-phonon 

coupled  system  and  admits  both  kink  soliton  and  polaron  solutions. 

Analytic  and  numerical  studies  show  that,  due  to  the  finite  response 

time  of  the  lattice  to  changes  in  its  dimerization,  the  kink  has  a 

maximum  free  propagation  velocity.   In  a  numerical  simulation  of  a 

simple  photoexcitation  experiment  in  trans-polyacetylene  we  find  that  a 

kink-antikink  pair  and  a  coherent  optical  phonon  package,  a  breather, 

are  produced.  We  calculate  the  optical  absorption  spectrum  of  the 

photoexcited  system  and  suggest  that  the  breather  may  account  for  the 

anomalous  subband  edge  absorption  in  trans-polyacetylene.  We  compare 

two  suggested  Hamiltonians  for  cis-polyacetylene  and  compare  their 


VI 


photoexcitation  dynamics  with  those  of  trans-polyacetylene.  We  further 
investigate  the  role  of  coherent  anharmonicity  in  simple  models  of  a 
finite  polyene  and  of  polyyne.   As  the  observed  transport  properties  of 
polyacetylene  depend  critically  on  the  presence  of  external  impurities 
we  study  the  dynamics  in  the  presence  of  model  defects  and  we  find  that 
a  single  kink  may  be  trapped  at  a  single  site  or  confined  to  a  short 
segment  of  the  system.   In  photoexcitation  experiments  on  the  defected 
system  we  find  that  kinks,  breathers,  polarons,  excitons  and  trapped 
kinks  may  be  produced. 
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CHAPTER  I 
INTRODUCTIOM 


Polymers  have  long  been  valjed  as  plastics  for  their  low  cost  and 
lightweight,  durable,  flexible  structure.  The  observation  in  1977  that 
by  adding  small  quantities  of  metallic  and  non-metallic  dopants  the 
conductivity  of  polyacetylene  can  be  varied  over  twelve  orders  of 
magnitude,  from  that  of  a  good  insulator  to  that  of  a  fair  metal,  has 
not  suprisingly,  produced  a  great  deal  of  interest.  Most  of  the  ensuing 
activity  has  been  in  four  broad  areas-:   technological  applications, 
synthesis,  physical  testing  and  characterization,  and  theoretical 
modeling. 

Amongst  the  suggested  technological  applications  for  polyacetylene 
and  related  materials  are  Schottky  diodes,  solar  cells,  lightweight 
metals,  and  high  power  storage  batteries. 

Ito  et  al.-^  first  synthesized  self-supporting  films  of 
polyacetylene  by  the  direct  polymerization  of  acetylene  gas  in  the 
presence  of  catalysts.   This  method  is  essentially  uncontrolled  and 
produces  samples  of  widely  varying  structure  and  physical  properties. 
Since  then  much  effort  has  been  made  to  produce  samples  with 
predictable  and  reproducible  physical  characteristics,  as  well  as  the 
synthesis  of  forms  of  polyacetylene  that  are  air  stable  and  soluble  in 
organic  solvents.   Another  area  of  major  interest  has  been  in  the 
synthesis  of  other  conducting  polymers.   Some  conjugated  polyaromatic 
heterocyclic  polymers  and  some  conjugated  block  copolymers  can  be 
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synthesized  using  conventional  organic  techniques.  Their  conductivities 
are  only,  however,  at  most  1-1 0!o  that  of  polyacetylene. 

The  full  barrage  of  standard  solid  state  spectroscopies  has  been 
used  to  characterize  polyacetylene":   electron  and  X-ray  diffraction 
have  been  used  to  study  its  structure;  Raman  scattering  to  determine  its 
conjugation  length;  infrared  and  optical  studies  of  its  excitation 
spectrum;  conductivity  measurements  of  its  charge  transport  properties; 
and  E3R  and  susceptibility  measurements  of  its  magnetic  properties.   The 
extreme  sensitivity  of  the  morphology,  structure  and  transport 
properties  of  polyacetylene  as  currently  synthesized  has,  to  some 
extent,  been  responsible  for  the  widely  differing  experimental  data  and 
interpretation.   However  one  particularly  significant  conclusion  has 
emerged:   ESR  experiments  have  shown  that  the  charge  carriers  in 
polyacetylene  are  spinless  and  therefore  are  not  electrons. 

The  theoretical  work  has  been  in  three  general  areas.  Ab  initio 
calculations,^  usually  at  the  Hartree-Fock  level,  give  a  good 
understanding  of  the  structure  and  properties  of  the  groundstate.  They 
cannot,  however,  even  qualitatively  describe  the  properties  of  the 
excited  system.  Phenomenological  Hamiltonians    that  include  strong 
electron-electron  interactions  can  well  describe  both  the  groundstate 
and  the  excited  state  system.  Analytic  calculations,  however,  are 
extremely  difficult  and  computational  difficulties  limit  studies  to  the 
static  properties  of  short  chain  systems.   In  1979  Su,  Schrieffer  and 
Heeger  (SSH)^~   proposed  a  simple  tight-binding  electron-phonon  coupled 
model,  which  naturally  explained  the  reverse  spin-charge  relation  as 
arising  from  kink  solitons.  This  static  discrete  model  and  its 
continuum  limit   *  -^  proved  to  be  analytically  tractable  and  their 


properties  have  been  extensively  investigated.   The  success  of  this 
tight-binding  model  is  particularly  remarkable  as  neither  electron- 
electron  interactions  nor  quantum  lattice  effects  are  included.   Each  of 
these  is  separately  difficult  to  include  and  studies  continue. 

None  of  the  above  theoretical  viewpoints  has  included  dynamics, 
which  are  clearly  essential  if  the  unusual  transport  properties  of  the 
conjugated  polymers  are  to  be  understood.   Here  we  shall  undertake  a 
systematic  study  of  the  dynamics  of  the  nonlinear  excitations  in  the  SSH 
model.   The  focus  will  be  in  two  areas:   first,  to  obtain  a  basic 
understanding  of  the  dynamics  of  nonlinear  excitations;  and  second,  to 
simulate  real  experiments  and  make  qualitative  and  quantitative 
comparisons  with  data.   The  major  task  is  to  establish  whether  the 
simple  electron-phonon  model  can  even  qualitatively  describe  the 
dynamics  of  the  real  system.   In  particular  it  is  necessary  to  identify 
physically  plausible  production  mechanisms  for  nonlinear  excitations;  to 
establish  their  robustness  in  the  presence  of  deviations  from  ideality 
e.g.  damping  and  disorder;  and  to  investigate  their  transport  and 
thermodynamic  properties. 

In  Chapter  II  the  relevant  parts  of  the  theory  of  the  electron- 
phonon  model  are  outlined.  The  idea  of  nonlinear  excitations  is 
introduced  and  the  limitations  of  the  model  discussed.   In  Chapter  III 
the  numerical  algorithm  for  the  solution  of  the  discrete  SSH  model  is 
discussed.   It  is  shown  that  the  static  excitations  found  in  the 
continuum  analysis  are  good  approximate  solutions  to  the  discrete 
equations  and  that  their  optical  properties  can  be  understood  from  the 
continuum  calculations.  Analytic  and  numerical  studies  of  the  dynamics 
of  a  free  kink  show  that  the  kink  propagation  velocity  is  limited  by  the 


finite  response  time  of  the  lattice.   In  Chapter  IV  a  simple 
photoexcitation  experiment  in  trans-polyacetylene  is  modeleh.   It  is 
found  that  both  kinks  and  a  coherent  optical  phonon  package,  a 
"breather",  are  produced  in  such  an  experiment  and  that  this 
photoexcitation  scenario  is  consistent  with  the  experimental  data.   In 
Chapter  V  cis-polyacetylene  and  the  finite  polyenes  are  discussed.   In 
Chapter  VI  it  is  shov/n  that  a  wide  range  of  nonlinear  excitations  can  be 
produced  during  the  photoexcitation  of  disordered  systems.  Their 
robustness  is  established  and  their  transport  properties  discussed. 
Chapter  VII  introduces  the  polyynes  and  suggests  that  nonlinear 
excitations  may  be  important  here,  also.  Chapter  VIII  presents  our 
conclusions. 


CHAPTER  II 
BASIC  THEORY 


Trans- Pol yacetylene 

Polyacetylene  is  a  linear  chain  polymer  consisting  of  a  spine  of 
carbon  atoms  with  a  single  hydrogen  atom  bonded  to  each.       It  exists 
in  two  isomeric  forms:  cis-polyacetylene  and  the  energetically  more 
stable  trans-polyacetylene.   The  present  discussion  will  concentrate  on 
trans-polyacetyene .   (A  full  discussion  of  cis-polyacetylene  will  be 
undertaken  in  Chapter  V.)  Of  the  four  carbon  valence  electrons,  three 
are  in  a  bonds  to  a  hydrogen  atom  and  to  the  two  neighboring  carbon 
atoms.  The  remaining  electron  is  in  the  it  band,  which  is  half  filled. 
One  thus  expects  trans-polyacetylene  to  be  a  conductor.  However, 
undoped  trans-polyacetylene  is  a  semiconductor  with  a  bandgap 
2A  -  1.'4-1.6  eV.   This  gap  is,  at  least  in  part,  due  to  the  Peierls 
effect  ':  for  a  one  dimensional  system  any  non-zero  electron-phonon 
coupling  lowers  the  energy  by  inducing  a  gap  at  the  Fermi  level.   For  a 
half  filled  band  the  band  gap  at  wavevector  k  =  ir/Za  increases  the  unit 
cell  size  from  one  (CH)  unit  to  two,  thus  dimerizing  the  lattice:  i.e. 
there  is  an  alternation  of  long  ("single")  and  short  ("double")  bonds. 

This  symmetry  breaking  can  occur  with  a  single  bond  either  to  the 
right  (A  phase)  or  to  the  left  (B  phase)  of  an  even  site  carbon  (fig. 
2.1).  Thus  there  are  two  structurally  distinct  though  energetically 
degenerate  groundstates  in  the  long  chain  limit.   A  kink  soliton  can 
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interpolate  between  them.   Physically  this  kink  comes  from  the  joining 
of  a  section  of  A  phase  to  a  section  of  B  phase  producing  a  pair  of 
(say)  single  bonds  next  to  each  other  (fig.  2.2)  and  thus  a  localized 
electron  at  this  defect.  This  picture  of  a  kink  localized  over  a  single 
lattice  spacing  is,  however,  over-simplified;  rather  it  is  extended  over 
a  few  lattice  spacings  due  the  competition  of  two  effects:   on  the  one 
hand  the  interpolating  region  between  the  two  groundstates  has  a  higher 
energy  than  the  dimerized  state  and  thus  favors  a  narrow  kink;  on  the 
other  hand  the  stiffness  of  the  electron  gas  favors  a  wide  kink.   For 
"realistic"  polyacetylene  parameters  the  width  of  the  kink  is  2E,  ~    10- 
20a  (a  is  the  lattice  constant).  As  we  shall  see  this  large  width 
justifies  the  use  of  a  continuum  approximation. 

That  the  kink  has  an  associated  localized  electronic  level  can  be 
understood  from  the  following  argument.  Consider  a  neutral  chain  of 
2n+l  atoms  of  a  single  phase.  There  are  2n  bonds  joining  the  2n+1  (CH') 
units.  On  joining  the  two  ends  of  the  chain  to  form  a  ring  an 
additional  bond  is  added,  giving  2n+1  bonds.  Thus,  at  some  point  on  the 
chain  there  must  be  two  single  (say)  bonds  adjacent.   A  ring  with  an  odd 
number  of  sites,  therefore,  has  a  kink  in  it.  With  2n+l  sites  there  are 
2n+1  associated  energy  levels.  So  for  any  Hamiltonian  that  has  charge 
conjugation  symmetry  (i.e.  for  every  electronic  level  at  energy  e  there 
is  a  level  at  energy  -e)  there  is  a  level  exactly  at  midgap  associated 
with  the  kink. 

There  are  2n  electrons  in  the  valence  band  that  pair  to  give  zero 
spin.  For  a  neutral  kink  the  midgap  state  is  singly  occupied  and  thus 
the  kink  has  spin  1/2.  Adding  an  electron  (hole)  gives  a  kink  charge  of 
-1  (+1)  and  spin  0.  These  novel  spin-charge  relations  provide  a  natural 


explanation  for  the  simultaneous  observation  of  high  electrical 
conductivity  and  low  Curie  susceptibility  :  the  charge  carriers  in 
doped  polyacetylene  are  charged,  zero  spin,  kinks. 


The  SSH  Model 
The  SSH  Hamiltonian 

Su,  Schrieffer  and  Heeger^  (SSH)  have  proposed  a  simple  Hamiltonian 
for  trans-polyacetylene  based  on  a  number  of  reasonable  assumptions: 

1 .  All  many-body  effects  can  be  incorporated  in  a  single  particle 
Hamiltonian. 

2.  The  effects  of  a  electrons  can  be  accounted  for  in  the  chain 
cohesion  and  lattice  dynamics.  X-ray   and  Nf'lR  '^  data  show  that 
a  -  1.22A  and  Ug  -  0.03±  0.01  A  (a  is  the  lattice  constant  and  Uq 
is  the  dimerization) .  Thus  excursions  of  a  (CH)  unit  from  its 
Bravais  lattice  site  are  small  and  the  harmonic  approximation  is  a 
useful  approach  to  the  lattice  dynamics.   The  ir  electrons  are 
highly  mobile  and  delocalized.  They  must,  therefore,  be  treated 
explicitly. 

3.  The  atomic  displacements  can  be  treated  as  classical  variables 
and  electron  states  can  be  evaluated  with  respect  to  fixed  atomic 
displacements.   The  adiabatic  approximation  can  be  expected  to  be 
valid  if  the  ratio  of  the  sound  velocity  to  the  Fermi  velocity  is 

small.   For  "realistic"  polyacetylene  parameters  this  ratio 

on 
is  ~  0.05.   Brazovskii  and  Dzyaloshinskii   have  shown  rigorously 

that  the  adiabauic  approximation  requires  y  <<  1 , 

where  p  =  Kw  /A  .   For  trans-polyacetylene  \i  ~    0.1. 


H.     A  linear  combination  of  atomic  orbitals  can  be  used  as  a  basis 
set. 

5.  Only  matrix  elements  between  states  with  bonded  sites  need  be 
considered  (tight  binding  approximation). 

6.  The  system  is  one-dimensional  and  interchain  effects  can  be 

1  Q 

neglected.   Experiments  suggest  that  the  interchain  overlap 

is  ~  100  times  less  than  the  intrachain  overlap  in  well  oriented 

samples.   The  effects  of  Interchain  coupling  and  cross-linking  may, 

however,  be  important  in  many  materials. 

The  SSH  Hamiltonian  is 

n  n 

where  M  is  the  mass  of  a  (CH)  unit,  K  is  the  force  constant,  t^  is  the 
intrinsic  electron  hopping  matrix  element,  a  is  the  electron-phonon 
coupling  constant.   C  (C  )  creates  (annihilates)  an  electron  on  site  n; 
Uj^  is  the  displacement  of  the  nth  (CH)  unit  ^rom  its  lattice  site.   The 
first  term  is  the  lattice  kinetic  energy;  the  second  is  the  lattice 
elastic  energy  (due  to  a  electrons);  the  third  is  the  -tt  electronic 
energy,  consisting  of  intrinsic  and  phonon-assisted  contributions. 


q 
Dimerized  Groundstate^ 

The  SSH  Hamiltonian  has  two  degenerate  uniform  solutions 


u^  =  ±(-1)"  u^.  (2.2) 


Now  defining  a  "staggered  order  parameter", 


u^=  (-l)"V'^o  ^2.3) 


then 


u  =  ±1  (2.4) 

o 


It  is  interesting  to  note  that  .vith  two  degenerate  groundstates  and  with 
a  local  energy  maximum  at  u=0  the  SSH  model  is  topologically  identical 
to  the  $  -model.  This  similarity  will  be  exploited  at  various  times, 
but  it  must  be  remembered  that  $   is  a  lattice  theory  whilst  in  the  SSH 
model  the  electrons  play  an  important  role. 

Certain  analytic  results  have  been  obtained  for  the  SSH  model.   In 
particular  it  can  be  shown  in  the  weak  coupling  limit  that  the  band 
gap,  2A^,  is  given  by 


2A^  =  8au^.  (2.5) 

o      o 


The  major  problem  of  the  SSH  model  is  that  very  few  exact  results 
can  be  obtained.   In  particular  no  analytic  form  for  the  kink  soliton 
has  been  found.   Numerical  studies,  however,  show  that  for  the 
variational  kink  profile,  A(x)  =  A  tanh(x/F)  (the  exact  0   kink 
solution),  the  kink  half  width,  f,  is  -  7-lOa.   This  suggests  that  a 
continuum  approximation  to  the  SSH  Hamiltonian  may  provide  a  useful 
basis  for  further  analytic  studies. 
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The  TLM  Model 
The  TLM  Haiiiltonian 

The  continuum  approxijiation  is  valid  if  the  electron-phonon 

2 
coupling  is  weak:   ,\  <<1  ,  where  X  =  2a  /irt  k.   For  SSH  parameters 

X  -  0.2.  The  fermion  operators  C  and  C  can  be  considered  as  the  sum 

n      n 


1  p  1  o 

of  left-  and  right-going  waves .' ^ 


n  1  r   /  N   ik„na   .   ,  ,   -ik^nai  ,„  ,. 

C   =  —  |u(n)e  F   -iv(n)e   F    .  (2.6^ 


ns     /—     S  S 

/a 


For  the  half-filled  band  the  system  dimerizes  to  give  a  lattice  period 
2a.   Thus  we  define 


[+,    \   2ik„na   ,,  ,   -2ik„naT  .„  „> 

u  =  -r—  [A  (n)  e   F   ^-  A(n)  e    F   J.  (2.7) 

n   Ua 


To  order  a,  and  dropping  the  rapidly  varying  terms  of  the  type  (-1)'^  the 
SSH  Hamiltonian  becomes  the  classical  TLM  (Takayama,  Lin-Liu,  Maki) 
Hamiltonian. 


2 


-A5  fdxA^(x)  +  Jdx  ¥*(x)  [-iv^o^  I-  +  A(x)oJy(x)    (2.8) 
^  2  •"  •'  F  3  oX        1 


where 


,„,  .    (-u(x)^    2   4K      .,  cai1/2       „.         -_  _, 
nx)    =  (^^^^]  ,  0)^  =  — ,  g  =  ta(j^J    .  Vp  =  2t^a       (2.9) 


and  0.  is  the  ith  Pauli  matrix. 
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Varying  eq.   2.8  with  respect  to  u  (x)  and  v  (x)  gives 


oU(X)     ./\/\         /     \ 

-IV^  -t: +  A(x)v(x)  =  £  u(x) 

F   dX  V 


+  iv^  ^^   +  A(x)u(x^  =  e  v(x) 

F   dX  V 


(2.10) 


These  are  a  pair  of  coupled  Dirac  equations  for  massless  particles  with 
eigenvalues  that  are  unbounded  below.  As  the  Fermi  sea  of  electrons  is 
not  appropriate  for  modeling  a  bounded  tt  band  an  artificial  energy 
cutoff  W  is  introduced  to  simulate  the  band  edges.  This  is  expected 

only  to  alter  the  energy  scale. 

+ 
Varying  eq.  2.8  with  respect  to  A  gives  the  consistency  equation 


A(x)  =  -g  /wg   Z'  [u  (x)  v(x)  +  V  (x)  u(x)]        (2.11) 

V  ,3 


v;here  the  sum  is  over  spin  components  and  occupied  electron  levels. 
Defining  the  new  functions 

f^.(x)  =  u(x)  +  iv(x)  f_(x)  =  u(x)  -  iv(x)     (2.12) 

equations  2.10  become 


3f_ 

£  f  (x)  =  -iv^  — i  (x)  ±  iA(x)  f_(x)  (2.13) 

V  ±         F  3x  + 


and  then 


2 

[2  8'"  2      9A(x)     2,  n1  ^  /  N    „         / -,  1  I,  s 

v„  — -   +  e   ±  V  -r- Ax)Jf(x)  =  0.        (2.14) 

F   „  2  \)    F   dx          ■'   ± 
9x 
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Thus  the  problem  has  reduced  to  the  solution  of  a  "Schrndinger-li'Ke" 
equation,  with  the  potential  being  a  function  of  the  order 
parameter  A(x).   Of  course  the  consistency  conditions  (2.11)  must  also 
be  obeyed.  Solutions  are  particularly  easy  to  obtain  for  the  class  of 
potentials  with 


''^^^    ^-^fT.-'I-  ^2.15) 


Two  particular  solutions  that  fall  into  this  class  are  the  purely 
dimerized  lattice  and  the  kink  solution. 


The  Dime r ized  Lattice" 

For  the  purely  dimerized  lattice:   A(x)  =  A  one  gets 


2 

t^F^^  "  %^  -  0    '.^^^    =  0-  ^2.16) 

3x 


This  has  solutions 


. ,         -v„k   +   CO  +  A 
/    N        ^■        ikx    r      F  O-i 

u    (x)    =  N,     e  I ; 

n  k  '-  v„k  ■' 


-V   k   -  oj  -  A 
v^(x)   =  N^  e^"^   [~-7-k -^  ^2.17) 


With  'k  =  2  7=  ^-;r-^     ' 

/2tt 


These  plane  wave  eigenfunctions  have  the  energy  dispersion  relation 
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2    2    ^2 

.u)   =  A   +  v-k  (2.18^ 

O     r 


where  the  consistency  equation  is  satisfied  by 


Aq  =  W  e   "\  (2.19) 


Here  one  sees  that  W  merely  sets  the  scale  of  energies, 


The  Kink  Soli  ton^  ^ '^'' 

There  is  also  a  kink  solution  (fig.  2.3a)  with  order  parameter 


A(x)  =  ±A  tanh(x/5)  (2.20) 

o 


where  f   =   v^/A  .   (The  minus  sign  gives  the  antikink  solution.)  This 
^    F  o  o  u 

has  both  plane-wave  solutions  and  a  single  localized  level 

at  e  =0  with  wavef unctions 
o 


u  (x)  =  N  sech(x/5) 


v  (x)  =  -iN  sech(x/5)  (2.21) 


and  N  =  [a  /^l  v^]^^^. 

o     o    F 


This  electronic  level  is  localized  on  the  kink  with  the  same  coherence 
length  as  the  lattice  distortion. 
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The  kink  creation  energy   '   can  also  be  caiculated  as 


E,  =  2A  /tt.  (2.22) 

k     o 


Although  the  TLM  model  is  neither  Galilean  nor  Lorentzian  invariant,  the 
kink  mass  can  be  calculated  in  the  ansatz  x(t)=x-Vt  by  equating  the  mass 
to  the  coefficient  of  the  (1/2)V  term.   This  gives  the  mass  of  the 
kink,  M,^  (where  M  is  the  mass  of  a  (CH)  unit) 

M,.  =  (|](^]'m  (2.23) 

for  SSH  parameters 


M  ~  3-6m  (2.2n) 


where  m^  is  the  electron  mass. 

Thus  the  kink  is  very  light  and  kink  dynamics  can  be  expected  to  be 
important. 


Th 


e  Polaron^^ '^^ 


It  is  natural  to  look  for  a  static  solution  to  the  TLM  equations 
that  does  not  involve  a  change  of  phase.  The  polaron  is  formed  by  the 
self-trapping  of  an  added  electron  or  hole  into  the  dimerized  lattice 
(fig.  2.3b).  The  polaron  is  topologically  equivalent  to  a  KK  pair,  as 
is  clear  from  its  order  parameter 


A(x)  =  A  -  k  v„  [tanh  k  (x+x  )  -  tanh  k  (x-x  ^1 
o    o  F  ■-      o    o         o    o  -' 
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2    2  2    2 

where  uj  =  k  v„  +  A  (2.25'> 

o    o  F    o 


and  tanh  k  x  =  (A  -  w  )/k  v„ 

o  o     o    o   Or 


This  self-consistency  condition  is,  however,  only  satisfied  in  trans- 
polyacetylene  if 


k  v^  =  A  //2.  (2.26) 

o  F    o 


This  lattice  deformation  is  supported  by  a  pair  of  intragap  levels 

symmetrically  about  the  Fermi  level,  which  in  trans-polyacetylene  are 

at  u)  =  ±A  //2.  For  an  electron  polaron  the  lower  level  is  doubly 
o     o 

occupied  and  the  upper  is  singly  occupied.  For  the  hole  polaron  the 
lower  level  is  singly  occupied  and  the  upper  is  empty. 
The  polaron  creation  energy  is 


E|^  =  2/2  Aq/tt  (2.27) 


i.e.  the  polaron  rest  energy  is  greater  than  that  of  a  single  kink,  but 
less  than  that  of  two  kinks  or  the  energy  (a  )  of  the  added  electron 
(hole). 


Limitations  of  the  SSH  and  the  TLM  Models 

Although  the  SSH  and  TLM  models  may  be  good  zeroth  approximations 
to  the  real  behavior  of  polyacetylene,  a  number  of  potentially  important 
effects  are  not  explicitly  included. 
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Electron-Electron  Interactions 

Experiments  shov;  the   first  excited  state  of  the  finite  polyenes 
((CH)^,  n=2,3,.,.)  has  ^A„  symmetry.   Band  structure  calculations 
predict  it  to  be  B  .  This  qualitative  disagreement,  and  its 
explanation  as  a  breakdown  in  the  single  particle  model,  has  been 
attributed  to  the  effects  of  electron-electron  interactions.   By 
analogy,  in  polyacetylene  the  single  particle  model  cannot  be  expected 
to  be  valid.  On  the  other  hand  certain  other  properties  of 
polyacetylene  are  well  explained  by  a  single  particle  model,  e.g.  the 
energies  of  the  infra-red  absorption  peaks  are  correctly  predicted."^ 

A  number  of  studies  on  electron-electron  effects  in  a  Peierls 
distorted  phase  have  been  undertaken  with  conflicting  conclu- 
sions.'-'" '       Calculations  at  the  Hartree-Fock  level  generally  show 
that  the  dimerization  is  decreased  by  electron-electron  effects,  whilst 
calculations  that  go  beyond  Hartree-Fock  conclude  that  the  dimerization 
is  increased. 

Although  electron-electron  effects  are  not  considered  explicitly, 
they  are  included  in  an  average  way.   The  fermion  operators,  as  in 
Landau  ferrai  liquid  theory,  can  be  considered  as  "fermion  quasi- 
particle"  rather  than  simple  electron  operators.  Also  the  values  of  the 
various  parameters  of  the  SSH  model  are  deduced  from  experimental  data 
and  thus  can  be  considered  as  being  renormalized  to  include  the 
electron-electron  effects.  However,  without  four-body  operators  in  the 
Hamiltonian  some  electron-electron  effects  are  evidently  omitted. 

Given  the  renormalization  already  considered  and  with  Horovitz's 
demonstration  that  the  optical  properties  of  the  SSH  model  are 
unaffected  by  the  addition  of  electron-electron  effects  (see  Chapter 
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III)  the  assumption  that  the  nonlinear  dynamics  and  statistical 
thermodynamics  of  polyacetylene  can  be  considered  as  mainly  arising  from 
the  electron-phonon  coupling  seems  to  be  justified.   All  many-body 
effects  can,  therefore,  be  considered  as  perturbations  to  the  single 
particle  Hamiltonian. 

Largely  because  of  the  controversy  over  how  to  include  electron- 
electron  effects  no  single  parameterization  of  the  SSH  model  is  accepted 
by  all  authors.^  Here  two  parameterizations  are  chosen.   "SSH 
parameters"^  are  chosen  to  reproduce  three  experimental  observables  of 
trans-polyacetylene:  the  bandgap  (2  A  ~  1.4eV),  the  dimerization 
(Uq  ~  0.014A)  and  the  full  bandwidth  (W  -  lOeV).   The  spring  constant 
(K  ~  21eV/A'")  is  the  value  for  ethane.   This  uniquely  fixes  the  electron 
phonon  coupling  constant  (ct  ~  4.1eV/A)  and  gives  the  coherence  length 
5  ~  7a.  For  many  of  the  simulations  performed  the  parameterization  W  = 
lOeV,  2A   =  3.8eV,  u^  =  0.1A,  C  =  2.7a  will  be  used.   These  "rescaled 
parameters"  "  are  numerically  more  convenient  as  they  enable  the 
consideration  of  systems  with  fewer  sites.   Also,  where  data  for  both 
sets  of  parameters  are  compared,  the  effects  of  discreteness  and 
acoustic  phonons  can  be  studied.   It  will  be  shown  that  there  are  no 
qualitatively  important  differences  in  the  numerical  results  between  the 
two  parameter  sets.   Of  the  quantitative  differences  none  is  important 
on  experimentally  relevant  timescales. 

Quantum  effects 

The  lattice  is  treated  classically.  The  two  phases  are  separated 
by  a  potential  barrier  and  are  classically  stable.   However,  quantum 
mechanically  their  wavef unctions  overlap  and  there  is  a  finite 


transition  rate  between  them.   Explicit  calculations  show  that  the 
dimerization  is  reduced  by  -  15-25?.3"»J^   Quantun  fluctuations  also 
reduce  the  kink  creation  energy  by  -  25?.-^   One  can  consider  these 
effects  as  being  included  in  the  renormalized  SSH  parameters.  Hov;ever 
there  are  certain  effects  that  can  only  be  understood  from  a  full 
quantum  treatment.   For  example,  classically  a  kink  may  be  pinned  to  a 
discrete  lattice.  This  is  not,  however,  experimentally  observable,  even 
in  principle,  because  of  the  quantum  uncertainty  in  the  kink  location. 
Similarly,  non-adiabatic  effects  may  have  important  consequences  for 
electronic  transitions,  absorption  coefficients,  phonon  absorption  and 
transport  properties. 

Interchain  coupling  and  intrinsic  disorder 

The  model  treats  each  polyacetylene  chain  as  independent.   However, 
X-ray  data   show  that  there  is  three  dimensional  ordering  in  trans- 
polyacetylene  with  the  antiparallel  alignment  of  dimerization  patterns 
on  neighboring  chains.   A  kink  introduces  a  defect^^  in  this  alignment 
with  an  energy  that  increases  linearly  with  the  separation  of  kinks. 
Thus  one  expects  a  KK  pair  to  be  confined  (though  not  necessarily  on  the 
same  chain)  to  reduce  the  misalignment  energy.  Experiments  suggest  that 
this  confinement  energy  is  -  3-IOK  and  that  a  KK  pair  is  confined 
over  ~  TOO  lattice  sites.   (As  will  be  seen  in  Chapter  V  this 
confinement  is  similar  to  the  intrinsic  confinement  in  cis- 
polyacetylene.  )  Interchain  hopping-^  of  charge  between  kinks  has  also 
been  suggested  as  the  dominant  conduction  mechanism  for  electrical 
conductivity  at  low  dopant  concentrations. 
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Other  morphological  effects  that  have  been  omitted  are  partial 
crystalinity ,  incomplete  isomerization ,  variations  in  the  chain  length, 
cross  linking  and  the  proximity  of  chain  ends  to  a  neighboring  carbon 
chain. 

Finite  Temperature  Effects 

Both  the  SSH  and  the  TLM  models  are  zero  temperature  models.   As 
the  band  gap  of  trans-polyacetylene  is  l.MeV  one  expects  thermal 
production  of  electrons  into  the  conduction  band  to  be  negligble. 
However,  the  self-focussing  of  thermal  energy  into  nonlinear  phonon 
wavepackets ,  breathers,  may  be  important. 

Extrinsic  Doping 

The  most  striking  experimental  feature  of  trans-polacetylene  is 
that  its  conductivity  can  be  varied  over  thirteen  orders  of  magnitude  by 
extrinsic  doping. ^-^   It  is  therefore  essential  that  the  system  be 
studied  in  the  presence  of  defects.  This  will  be  the  subject  of  Chapter 
VI. 

Dynamics 

If  the  mechanism  for  electrical  conduction  in  polyacetylene  is  to 
be  understood,  it  is  clearly  necessary  to  understand  the  dynamics  of  the 
nonlinear  excitations  of  the  model.  Attempts  to  find  time  dependent 
nonlinear  solutions  to  TLM-like  equations  have  proven  unsuccessful  (the 
TLM  Hamiltonian  is  neither  Lorentzian  nor  Gaililean  invariant — there  are 
two  velocity  scales,  the  Fermi  velocity  and  the  sound  velocity). 
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All  of  the  above  deficiences  of  the  basic  model  need  to  be 
addressed.  A  first  logical  step  is  to  study  the  dynamics  of  the  basic 
3SH  model.   This  is  the  subject  of  Chapter  III. 


Figure  2.1 
The  tv;o  degenerate  phases  of  trans-polyacetylene , 
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CHAPTER  III 
DYNAMICS  OF  A  SINGLE  KINK  IN  THE  SSH  AND  TLM  MODELS 


In  this  chapter  a  systematic  study  of  the  dynamics  of  a  single  kink 
in  both  the  TLM  and  SSH  models  will  be  pursued.   Before  this  is  done  it 
is  useful  to  outline  the  numerical  method  employed  and  discuss  boundary 
conditions.   A  brief  study  of  the  optical  absorption  of  some  simple 
static  lattice  configurations  shows  the  similarity  in  nonlinear  behavior 
of  the  SSH  and  TLM  models. 


on 
The  Numerical  Tecnnigue  ^ 


The  energy  of  the  SSH  Hamiltonian  can  be  formally  written  as 


"sSH  =  i^  %^'^nu^}]  (3.1) 

n 


where  the  first  term  is  the  lattice  kinetic  energy  and  the  second  is  the 
potential  energy  (V),  which  is  a  functional  of  the  atomic  displacements 
u^.  The  potential  energy  is 


V  =  •:r  Z  (u  -  u   J   +   I  m    e  (3.2) 

2n    n-1  v,sv,s 

n  v,s 


where  the  first  term  is  the  lattice  strain  energy  and  the  second  is  the 
lattice  electronic  energy.  The  sum  is  over  occupied  states  of 
energy  e  and  occupation  m.  The  electronic  energy  can  be  found  by  the 
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direct  diagonalization  of  the  hopping  matrix  T,  which  has  components 


^ij  =  l^o^"t^-'^j]M5,,j_T  ^^i,j.i}.  (3.3) 


This  has  eigenvalues  e     given  by 


(T  -  E  I)  =  0 


where  E  =  (e, ,e„  .  .  .  .  e„) 

—  \       d  N 


Newton's  equations  of  motion  are 


(3.4) 


5v[iu  n 

MQ  =  -  ——-^^—   .  (3.5) 


n      ou 
n 


The  time  derivative  can  be  evaluated  by  finite  differences  as 


%  (f")  =  -^t  [^n^'"^  "  -"n^""^^^  ^3.6) 


where  dt  is  the  time  between  the  m-1  and  the  m  iterations.  This  gives 
the  equation  of  motion  as 


2 
u^(ra)  =  2u^(m-1)  -  u^(m-2)  -  l^liL  |^(|  ^^(ni-i  ) }  ]  .      (3.7) 

n 


The  functional  derivative  of  the  potential  can  also  be  calculated  by 
finite  differences  as 
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6V[{uJ) 

r^ =  [v(u  u  .  .  .u  +5u. .  .u  )  -  V(u  u  .  .  .u  .  .  .u  )]/6u         (3.3) 


where 


6u  <<  u^.  (3.9) 


Energy  conservation  is  used  as  a  check  on  the  stability  of  the 
algorithm.   In  general  energy  is  conserved  to  better  than  99.99?. 

At  some  points  in  this  study  spatially  homogeneous  velocity 
dependent  damping  is  applied  by  the  addition  of  the  term 


-^^n^S^n  ^3.10) 


to  the  equations  of  motion.   This  gives  evolution  equations 


.  ,  ,2  6V[{u  }) 

^n  =  TTT-T  [(2+M)u^(m-1)  -  u^(m-2)  -  ^ TT^]  •    ^3.1^ 

n   (l+u)        n        n         M      6u 

n 


Also  in  some  studies  instantaneous  damping  is  applied  by  setting 
the  velocity  to  zero  at  every  timestep,  i.e., 


u  (m-1)  =  u  (m-2).  (3.12) 

n        n 


This  modifies  the  equations  of  motion,  eq.  3.7,  to 


u^(m)  =u^(m-1)-l^^!^.  (3.13) 

n 
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Boundary  Conditions 

Two  types  of  boundary  conditions  are  used  in  this  numerical 
study. 3S 

Ring  Boundary  Conditions 

A  ring  is  a  system  in  which  there  is  a  finite  coupling  between  the 
1st  and  the  Nth  sites  of  a  system  of  N  lattice  sites.   Ring  boundary 
conditions  are  imposed  by  setting 


%  -   Vn  ^3.14) 


and 


C  =  C,   .  (3.15) 

n    N+n 


As  discussed  in  Chapter  II  a  ring  with  an  odd  number  of  sites  is 
topologically  restricted  to  contain  a  single  kink  (or,  more  generally, 
an  odd  number  of  kinks).  If  the  kink  (antikink)  is  assigned  a 
topological  charge  of  +1  (-1)  then  an  odd  length  ring  must  have 
topological  charge  of  ±1 .  On  the  other  hand  a  ring  with  an  even  number 
of  sites  must  contain  an  even  number  of  solitons  (including  zero).   It 
thus  has  topological  charge  0.  This  difference  between  even  and  odd 
length  rings  persists  for  all  values  of  N  and  leads  to  the  boundary 
condition  on  the  staggered  parameter 


u  =  (-1)    u  „. 
n  n+N 
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There  is  also  a  more  subtle  distinction  between  systems  of  length 

Mn  and  4n+2.  For  an  undimerized  system  (u  =  0)  the  eigenstates  of  the 

o 

SSH  Hamiltonian  are  Bloch  states  with  energy 


E  =  2t  cos  k  a  (3.17) 

n     o     n 


where  the  wavevector ,kj^,  satisfies  the  quantization  condition 


k^  =  2^n/N  ^  <  n  <  |.  (3.18) 


Thus  for  N=4n  and  a  half  filled  band  the  Fermi  level  lies  at  E=0.   By 
the  Peierls  effect  the  system  then  spontaneously  dimerizes.  On  the 
other  hand  for  N=4n+2  the  Fermi  level  lies  at  the  center  of  a  gap  of 
width  4t  sin  (tt/N).  Thus  the  4n+2  length  ring  only  dimerizes  for 
sufficiently  large  values  of  N.   (The  critical  value  of  N  depends  on  the 
coupling  constant,  A.) 

This  difference  between  4n  and  4n+2  length  rings  decreases  as 
1/N.  Most  of  the  numerical  studies  performed  here   are  for  N=98  where 
the  difference  is  small.   (Finite  N  dependence  is,  however,  interesting 
and  will  be  studied  briefly  in  Chapter  V.) 

Chain  Boundary  Conditions 

On  a  chain  there  is  no  direct  coupling  between  the  Ist  and  Nth 
sites  of  a  system  of  N  sites.  This  results  in  a  breaking  of  the 
electron-hole  symmetry  to  order  1/N,  which  for  all  reasonable  chain 
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lengths  is  small.   Here  two  different  types  of  chain  boundary  conditions 

will  be  used. 

1.  In  "free  boundary  conditions"  no  constraint  is  applied  to  the 
end  sites.   To  conserve  the  chain  length,  however,  it  is  necessary 
to  add  to  the  Hamiltonian  the  "pressure  term" 


IT     1       N 


2.   In  "fixed  boundary  conditions"  each  of  the  end  sites  is  bound 
to  its  nearest  neighbor  by  the  constraint 

^1  =  '^2'       '^N  =  Vr  ^3.19) 

At  no  point  in  this  study  has  any  qualitative,  or  significant 
quantitative,  difference  been  found  between  the  dynamics  of  chains  with 
free  and  fixed  boundaries. 

Statics  of  the  SSH  Model 


As  a  preliminary  to  the  study  of  the  dynamics  of  the  SSH  model  it 
is  important  to  verify  that  the  static  analytic  solutions  to  the  TLM 
Hamiltonian  well  approximate  the  solutions  to  the  SSH  equations. 

As  an  initial  condition  we  use  the  t=0  lattice  configuration  and 
minimize  the  energy  of  the  system  by  varying  the  dimerization  u^  from 
its  continuum  value.   The  shape  of  the  initial  profile  remains  fixed 
however.   The  solution  to  the  SSH  Hamiltonian  is  then  found  by 
numerically  relaxing  around  this  initial  configuration  (see  eq,  3- 13). 
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This  allows  the  shape  of  the  nonlinear  excitation  to  change  by  the 
emission  of  phonons. 

It  has  been  shown  analytically  that  the  dimerized  lattice  is  an 
exact  solution  to  the  SSH  equations:   the  energetically  minimized 
dimerized  groundstate  produces  few  phonons  and  appears  to  be  dynamically 
staole  for  all  times. 

Although  no  analytic  form  for  the  kink  in  the  SSH  model  exists  it 
is  expected  that  the  TLM  kink  profile  should  be  a  good  approximation. 
On  relaxing  the  kink  profile  there  is  a  strong  emission  of  acoustic 
phonons,  which  appear  in  the  staggered  lattice  distortion  as  a  "saw- 
tooth" (fig.  3.1a).   The  effect  of  the  acoustic  phonons  can  most  easily 
be  seen  by  plotting  the  short  and  long  range  components  of  the  staggered 
order  parameter  separately.   The  actual  ion  displacement  at  site  n  is 


u  =  (-1 )^  u  +  n.6a  (3.20) 


where  5a  is  the  change  in  the  lattice  constant  and  u  is  the  staggered 
order  parameter  with  respect  to  the  local  lattice  constant.  Now 
def  ine^''' 


r     =  I   (-1)"  [2u  -  u   ,  -  u   1  (3.21) 

n   4       '•  n    n+1    n-V  ^j.   ^ 


^n  =U^%   ^    Vl  ^  'Vll-  (3.22) 


Then  r^  _  j^  and  s^^  -  n.5a.   In  these  variables  the  kink  profile  is  more 
easily  seen  (fig  3.1b).  Note  that  the  dip  in  s^  corresponds  to  a 
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contraction  of  the  lattice  constant  around  the  kink  (Tor  proof  see 
Appendix  C  of  ref.  15). 

Optical  Absorption  within  the  SSH  Model 

Before  discussing  the  optical  absorption  within  the  SSH  model  it  is 
useful  to  review  the  optical  absorption  in  the  TLM  model.   '^°  The 
optical  absorption  is  defined  within  the  linear  response  dipole 
approximation  as 


_  A  - 

ctCto)  =  TT  f  —  A   E   |<y_|o,|-f,>|   6(E,  -  E,  -  w)       (3.23) 


where  a-,  is  the  current  operator,  the  6-f unction  ensures  the 
conservation  of  energy  and  A  is  an  unimportant  constant.   Figure  3-2 
shows  the  available  electronic  transitions  for  the  dimerized  lattice, 
the  kink  and  the  polaron.  Using  the  analytic  form  of  the  wavef unctions 
(eq.  2.17)  the  optical  absorption  of  the  purely  dimerized  lattice  can 
easily  be  calculated  as 


ap(a))  =  k[2^^/i^][J  -   Hh^]    '^^.  (3.2U) 


This  is  divergent  at  the  band  edge  and  decreases  as  the  inverse  square 
of  the  frequency. 

Great  care  is  needed  in  calculating  the  optical  absorption  of  the 
single  kink.^^  '^^  Although  the  kink  is  a  local  excitation  it  results  in 
a  global  change  in  the  dimerization  of  the  system.   (To  one  side  of  the 
kink  the  system  is  in  the  A  phase,  on  the  other  side  it  is  in  the  B 
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phase.)  This  alters  the  phase  relations  in  the  valence  and  conduction 
bands,  causing  diagonal  transitions  (Ak=0)  to  be  explicitly  forbidden. 
For  the  midgap  transition 

0^    ^         tt'^    „    5    r    2         ,2-1-1/2  v,2r    TT       (2  2^-1/21  ,^    ^^, 

aj,    (o))    =  Tj     A  ^   [co      -  A^j  sech   [^^   [co      -   A^J  J         (3.25) 

o 
whilst  for  the  interband  absorotion 


a^^(aj)  =  a^{^)    [l  -  25/l].  (3.26) 


It  is  important  to  note  that  the  kink  removes  weight  essentially 
uniformly  in  frequencies  above  the  band  edge  because  the  alteration  in 
the  dimerization  over  a  semi-infinite  region  requires  a  large  number  of 

Fourier  modes  for  its  description. 

21 
The  polaron  absorption   has  similarly  been  calculated  and  consists 

of  three  contributions  plus  the  interband.  Transitions  between  the  two 

intragap  polaron  levels  give  ct  (co).  There  is  a  second  contribution, 

P 
2 
a  (o)),  arising  from  transitions  between  the  valence  band  and  the  lower 

polaron  level  (for  the  hole  polaron)  and  between  the  upper  polaron  level 

and  the  conduction  band  (for  the  electron  polaron).  Contributions 

between  the  valence  band  and  the  upper  polaron  level  and  between  the 

lower  polaron  level  and  the  conduction  band  give  the 

contribution  a  (;o).   The  relative  intensities  of  each  of  these  is  shown 
p 

in  fig.  3.3  as  a  function  of  the  location  of  the  intragap  polaron 

level.   (For  the  polaron  in  trans-polyacetylene  w  =  1//2.   For  cis- 

polyacetylene  see  Chapter  V.)  For  most  gap  locations  a  (;u)  is  much 

P 

■i 
stronger  than  a  (u)  even  though  it  involves  only  one  transition, 

P 
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3 
whilst  ct  (o))  involves  many  transitions.   This  suprising  intensity  ratio 

is  again  a  result  of  a  fortuitous  phase  relation  in  the  continu'um 

equations. 

This  continuum  theory  suffers  from  a  number  of  limitations;  the 
most  important  of  which,  is  that  only  a  few  simple  lattice 
configurations  can  be  studied.  Further  it  is  not  clear  whether  the 
phase  relations  discussed  above  will  hold  in  the  more  realistic  discrete 
system. 

A  theory  for  the  optical  absorption  in  the  SSH  model  has  been 
constructed  by  Horovitz.^   A  review  is  given  here    (full  details  are 
given  in  Appendix  C).  Starting  from  the  definition  of  the  charge 
density,  p(n)  =  eC  C  ,  and  relating  it  to  the  current  density  through 
the  continuity  equation,  he  found  the  current  operator  j(n)  to  be  given 
by 

j(n)    =  -ie[t      +  ci[u     -   u      Jl[c'"C      ,    -  c""    ,C    1.  (3-27) 

•^  '-    o  ^    n         n+r-"-    n  n+1  n+1    n-'  ^j      i  ' 

Then  in  the  linear  response  dipole  approximation 


a(a))  =  ^e^t^  \  E   |M,  ,|^6(e,  -  e,  -  w)  (3.28) 

o  u)N      '  I  ,2 '     2     1 


where 


M,^2  =  '  t^  ^  F"  ^%  ~  Vl^l^M^^^^^"^^^  ~  f/n+l)f*(n)]    (3.29) 


and  f  is  the  wavef unction  of  the  a  level. 

a 

He  further  showed  that  a(a))  obeys  a  sum  rule 
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\      a(a))du)  =  -  -|^  Tre^  <\{^^>  (3-30) 

o 


where  H   is  the  electronic  energy  of  the  configuration.   This  algorithm 

for  the  absorption  is  easily  implemented.   It  is  important  to  note  that 

all  these  results  remain  unchanged  when  electron-electron  effects  are 

added  at  the  Hubbard  level,  i.e.  H   =  Z  V   p(n)  p(m)  is  added  to  the 

ee       nm 

nm 

SSH  Hamiltonian.   This  is  because  [h   ,  p(n)]  =  0. 

Figure  3.^  shows  the  numerically  calculated  optical  absorption  of 

the  purely  dimerized  lattice.   The  5-function  in  eq.  3-28  has  been 

approximated  by  a  Lorentzian  of  width  G  (for  all  the  present  studies 

G  ~  A  /Id.   The  sum  rule  is  obeyed  to  better  than  99.92?. 
o 

Figure  3-5  shows  the  optical  absorption  of  two  widely  spaced  kinks 
on  a  ring  of  98  sites.  As  expected  there  is  a  strong  midgap  absorption 
in  addition  to  the  interband  continuum.  Here  the  sum  rule  is  obeyed  to 
>  99.97?  .  In  fig.  3.6  the  difference,  Aa((o),  between  the  absorption  of 
the  ring  with  the  two  kinks  and  the  absorption  of  the  purely  dimerized 
ring  is  plotted.  Here  the  expected  interband  bleaching  is  evident  over 
all  frequencies  above  the  band  edge. 

Figure  3.7  shows  the  optical  absorption  of  a  single  polaron  on  a  98 

1         2 
site  ring.  The  contributions  a  ((d)  and  a  (o))  are  evident  at  --  3A  /? 

P        P  o 

and  -  A  /4  respectively.  The  weak  contrbution,  a  (to),  is  present  but 
o  p 

only  ~  2%   of  that  of  a    (oj).   Figure  3.8  shows  Aa(to)  for  the  polaron. 
Note  that  the  polaron  removes  intensity  over  a  wide  energy  range  above 
the  band  edge.  This  is  because,  for  these  parameters  (^  -  2.7),  the 
polaron  is  is  narrow  and  requires  a  large  number  of  Fourier  modes  for 
its  description. 
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These  simple  calcjlations  show  that  the  optical  absorption  of  these 
static  configurations  is  in  good  agreement  with  the  continuum 
calculations.   The  valiiity  of  the  sum  rule  is  also  confirmed.   In 
subsequent  chapters  this  formalism  will  be  used  to  calculate  the  optical 
absorption  spectrum  of  various  dynamical  excitations  of  the  system. 

Kink  Dynamics  in  the  TLM  Model 

It  was  noted  in  Chapter  II  that  the  TLM  equations  are  neither 
Galilean  nor  Lorentzian  invariant.   There  is  thus  no  obvious  velocity 
scale  in  this  problem.  However  due  to  the  finite  response  time  of  the 
lattice  to  changes  in  its  dimerization,  the  kink  has  a  maximum  free 
propagation  velocity  as  can  be  seen  from  the  following  simple 
argument  .^' '-^^ '    For  a  kink  of  velocity  v  the  lattice  changes  from  the 
A  to  B  phase  (say)  over  a  distance  ^(v)  and  in  time  ^(vVv.  The 

minimum  response  time  of  the  lattice  is  the  inverse  of  the  renormalized 

1/2 

phonon  frequency,  a)„,  where  u)„  =  (2X)   w   (w  being  the  bare  phonon 

H         K  o   o 

frequency).  This  gives  the  maximum  kink  propagation  velocity  (v  ) 

as  V  ~  aj„  Efv  1  .   Assuming  that  the  kink  width  is  not  strongly  velocity 
m    R    m' 

dependent,  then 


v^/v^  -  (2X)'^^  (5/a).  (3.31) 


For  SSH  parameters  [E,  ~   7a,  \  ~   0.2)  this  gives  (Vj^/v^)'"  -  19 — the 
maximum  velocity  is  greater  than  the  sound  velocity. 

The  conclusion  of  the  above  argument  can  be  confirmed  from  an 
approximate  analytic  calculation.   Using  the  ansatz  for  the  kink  profile 
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A(x,t)  =  A^  tanh  (x  -  vt)/^(v)  (3.32) 


it  is  straightforward  to  show  that  the  lattice  kinetic  energy  of  the 
kink  (T)  is 

(v/v  ^^ 
T  H  I  J  A  (x)dx  =  A^Y  ~-  (3.33) 


where 


,(v)/c(0)  (3. 34) 


and 


Y  =  A^/24TTAt^.  (3.35) 


It  has  been  shown  numerically^  that  the  soliton  width  is  velocity 
dependent  (contrary  to  the  assumption  of  the  simple  argument  above);  the 
excess  energy  due  to  the  lattice  deformation  being 


Ej^  =  8A^(1  -  r)^  (3.36) 


where  numerically  we  find 


,(5  =  7)  ~  0.20  (3.37) 

;(^  =  2.7)  -  0.10. 


Then  the  Langragian  of  a  kink  of  finite  velocity  is 
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L  =  A^lsd  -  rr   -  Yq^/r|.  (3.38) 


T'nence  the  momentum  (p  )  conjugate  to  the  displacement  q  (q  =  v/v  )  is 


Pq=^  (3.39) 


and  the  momentum,  p  ,  conjugate  to  the  width  r  is 


Pr.  =  0.  (S.'^O) 


r 


Now  extremizing  the  Lagrangian  with  respect  to  r  and  q  gives 


-Yq^  +  26r^(I  -  r)  -  0  {3.^^) 


'1-^-0.  (3.42) 


Then  for  consistency 


r  =  0,  q  =  0.  (3.^43) 

Thus  from  eq.  3.41  the  velocity  has  an  extremum  at  r=2/3  independent  of 
the  parameters  y   and  B. 

Substituting  back  for  r  the  kink  has  a  maximum  propagation  velocity 
(v^)  given  by 

[v/vj2  =  (|_)  g/Y  (3.44) 
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which  gives  the  momentjm  at  the  maximum  kink  velocity  as 


P^=|78.  (3.^5) 


Thus  the  energy  at  maximum  kink  velocity  is 


E^=|6A^.  (3.45) 


The  kink  energy  can  also  be  written 


^  ,ru..2. 


E|^(p)    =   A^    [p^/4Y   -   pV6itY^3].  (3.^7) 


This  gives  a  maximum  in  the  energy  at 


P^  =  8Y6.  (3.43) 


This  apparent  contradiction  between  eqs.  3.45  and  3.48  is  resolved  by 
noting  that  the  previous  calculation  was  for  the  momentum  at  the  maximum 
kink  velocity,  whilst  this  calculation  gives  the  momentum  at  which  the 
energy  is  a  maximum.  Now  defining  the  effective  kink  mass  by 


2 
1     fi  P 

dp 


gives 


if  =  k  -  16  P'/^S-  (3.50) 
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Thus  in  the  low  momentum  limit  the  kink  behaves  like  a  fJewtonian 
particle  of  mass  2  a.   However,  the  effective  mass  diverges  at  the 
critical  value  of  the  momentum 


2   8 
P(,  =  3  YB.  (3.51) 


This  is  the  same  as  the  kink  momentum  for  free  propagation  at  the 
maximum  velocity  (eq.  3.^5).  Using  appropriate  value  for  B  and  Y  one 
gets  for  S3H  parameters 


(v  /v  )^  ~  n,   m*  -  5.6  m  ,   E   -  0.1  A  (3.52) 

m  s  em       o 


and  for  rescaled  parameters 


(v  /v  )^  ~  1.2   m  ~  80  m  ,   E  ~  0.06  A  .        (3.53) 
m  s  em       o 


The  maximum  velocity  calculated  here  for  SSH  parameters  is  in  excellent 
agreement  with  that  calculated  above.  The  mass  of  the  kink  is  also  in 
good  agreement  with  that  calculated  in  Chapter  II, 

Kinks  Dynamics  in  the  SSH  Model^' 

To  test  the  validity  of  the  above  calculations,  a  detailed 
numerical  study  of  kink  dynamics  in  the  SSH  model  was  undertaken.   A 
single  kink,  given  by  its  continuum  analytic  form,  was  given  a  Galilean 
boost  towards  the  center  of  a  79  site  chain  with  fixed  boundary 
conditions.   (All  results  were  found  to  be  similar  for  free  boundary 
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conditions  and  for  chains  of  differing  lengths.)  In  fig  3.9  the  short 
time  kink  velocity  squared  is  plotted  as  a  function  of  the  energy 

input.   Both  axes  are  in  dimensionless  units,  where  unity  corresponds  to 

2 
the  theoretical  value  of  v  and  the  theoretically  calculated  energy,  E^, 

at  this  maximum  velocity.  The  solid  line  comes  from  the  continuum 

calculation  of  the  previous  section.   It  is  clear  that  the  data  for  the 

two  parameter  sets  lie  approximately  on  a  single  curve.   Indeed  by 

treating  B  and  Y  as  free  parameters  and  force  fitting,  better  agreement 

may  be  obtained.  At  low  energy  inputs  the  data  are  in  good  agreement 

with  the  continuum  theory.   The  maximum  Kink  velocity  is  also  well 

predicted  by  the  continuum  calculation.   However  the  continuum  theory 

strongly  underestimates  the  energy  input  required  to  reach  this 

velocity.  This  is  a  result  of  the  simplicity  of  the  ansatz,  which 

allows  for  kink  translation  and  uniform  contraction  only.  In  the 

numerical  studies  the  kink  is  also  coupled  to  the  optic  and  acoustic 

phonon  fields,  which  provide  efficient  routes  for  energy  dissipation. 

Indeed,  n-umerically  we  find  that  phonon  production  is  strongly  favored 

over  kink  contraction;  the  width  of  the  kink  only  decreasing  by  -10% 

even  at  the  highest  velocities. 

Short  time  dynamics  have  only  been  considered  up  to  this  point;  as 

after  approximately  one  phonon  period,  phonon  emission  becomes 

important.  Optic  phonons  couple  to  the  nonlinear  excitations  to 

order  (a/^)  whilst  acoustic  phonons  couple  to  order  (a/E,)    .   Thus  for 

rescaled  parameters  phonon  effects  should  be  much  stronger  than  for  the 

SSH  parameters.  Figure  3-9  also  shows  the  energy  input  as  a  function  of 

the  kink  velocity  squared  over  the  time  period  0.15-0.30  psecs ,  during 

which  the  kink  propagates  at  an  approximately  constant  velocity.  Here 
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the  p'nonon  effects  reduce  the  maximum  kink  velocity  by  a  factor  of  ~2, 
for  f~  7,  but  by  a  factor  of  -10  for  the  more  discrete  kink 
with  E.~   2.7.   Figures  3.10  (a)  and  (b)  compare  the  lattice  dynamics  of  a 
single  kink,  with  rescaled  parameters,  with  energy  inputs  S-0.3  E 
and  E~2.0  E  respectively.   At  low  energy  inputs  the  kink  propagates 
essentially  uniformly,  whilst  at  higher  energy  inputs  kinetic  energy  is 
lost  by  the  emission  of  an  optical  phonon  package.    This  may  be  viewed 
as  the  classical  analog  to  the  emission  of  Cerenkov  radiation  by  a 
particle  traveling  faster  than  its  maximum  propagation  velocity  in  the 
medium. 

Although,  as  we  shall  see  in  Chapter  VI  cnarge  conduction  by  the 
ballistic  transport  of  kinks  is  unlikely,  kink  propagation  is  essential 
to  the  understanding  of  the  short  time  dynamics  of  the  photoexcited 
system. 
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Figure   3- ' 

The  Single  Kink 

(a1  unsmoothed  order  parameter,  u   .   for  single  kink, 
(b)  Smoothed  order  parameter,  r  ,  for  the  same  single  kink  and  s^ 
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Figure  3.3 


Relative  strength  of  transitions  in  the  polaron  as  a  function 

of  the  location,  u  /A  ,  of  the  intragap  levels.  For 
,o  o  .  ,       , 


trans-polyacetylene  w  /A  -  0.7. 
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Figure   3'^ 

Numerically  calculated  optical  absorption,  a(w),  of  a  93  length 

dimerized  ring  with  rescaled  parameters.   The  full  band  gap 

2A   =  3.9  eV  and  the  full  bond  with  B  =  ^t^  =  10  eV.   The 

6-functions  in  energy  have  been  approximated  by  Lorentzian 

lineshapes  of  width  -  A  /lO.   In  this  and  subsequent  figures  the 

vertical  axis  is  in  arbitrary  (but  consistent")  units. 
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Figure  3.5 
Optical  absorption,  a(uj),  of  two  widely  spaced  kinks  on  a  98  site  ring. 
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Figure  3.6 

Difference  between  optical  absorption  of  two  kinks 
and  the  dimerized  lattice  of  a  98  site  ring. 
(Figure  3.6  =  Figure  3.5  -  Figure  3.^^). 
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Figure  3.7 
Ootlcal  absorption,  a(w),  of  a  polaron  on  a  98  :3ite  ring, 
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Figure  3.3 

Difference  in  absorption  of  a  single  polaron 
and  the  purely  dimerized  lattice  of  a  98  site  ring, 
(Figure  3.8  =  Figure  3.7  -  Figure  3.^) 
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CHAPTER  IV 
PHOTOEXCITATION  IN  TRAMS-POLYACETYLENE 


We  begin  this  study  of  photoexcitatlon  in  trans-polyacetylene  with 
a  simple  experiment:   '-^'  the  excitation  of   a  single  electron  from  the 
top  of  the  valence  band  to  the  bottom  of  the  conduction  band.  To 
simulate  an  infinite  system  98  sites  are  used  and  periodic  boundary 
conditions  applied.   In  the  rest  of  this  work  the  "rescaled  parameters" 
will  be  used  as  they  allow  smaller  systems  to  be  studied,  thus  reducing 
the  computation  time.   (The  dynamics  are  expected  to  be  similar  for  the 
more  realistic  SSH  parameters.)   At  t=0  a  single  electron  is  manually 
removed  from  the  top  of  the  valence  band  and  placed  in  the  bottom  of  the 
conduction  band.  All  electrons  then  remain  in  the  same  levels 
throughout  the  experiment  and  the  system  simply  evolves  adiabatically. 
As  electronic  transitions  should  only  be  important  on  the  nanosecond 
timescale  this  should  be  a  good  approximation  for  the  subpicosecond 
dynamics  considered  here. 

The  single  electron  photoexcitatlon  adds  energy  2A  to  the 

system.   In  less  than  0.1  psecs  a  kink-antikink  (KK)  pair  is  formed. 

The  kink  and  antikink  separate  with  their  maximum  free  propagation 

velocity,  v^,    (fig.  4.1).   Each  kink  has  creation  energy  2A  /tt  and 

kinetic  energy  <0.1  A  .   Thus  the  two  kinks  have  a  total  energy 

of  <1.5  A  .   Figure  4.1  shows  that  a  substantial  fraction  of  the 
-     o 

remaining  energy  from  the  photoexcitatlon  is  localized  as  a  nonlinear 
lattice  excitation — an  amplitude  breather.   The  breather  dynamics  were 
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studied  by  removing  the  central  42  atoms  from  the  photoexcitation 
experiment  at  t=0.301  psecs  and  allowing  them  to  evolve  separately  (fig. 
4.2a').   It  was  also  shown  numerically  that  the  breather  is  stable  with 
no  appreciable  change  in  frequency  or  amplitude  after  5  psecs. 

The  Analytic  Breather 

As  the  breather  is  formed  by  small  oscillations  about  the  purely 
dimerized  lattice,  one  expects  that  an  approximate  analytic  form  may  be 
obtained  from  the  known  effective  Lagrangian'^ '  for  the  half-filled  band 

'2    -2 

L  =  —  |l/2  A^[5Ln(2E  /A)  +  1/2]  -  ^^/^\    -   vJ   ^—  +  ^—  }     (14.1) 
■"■^TT  c         -■  F   T,,,2    „  2 

F  24a     2a)n 

a 

where  E^  is  the  electronic  cutoff  energy  and  the  dot  and  prime  denote 
time  and  spatial  derivatives  respectively.   (This  Lagranglan  cannot  be 
used  to  find  kink-like  solutions  as  it  is  ill-defined  as  A->'0.)  From  eq, 
4.1  the  equation  of  motion  is 


!  2  " 

A/2X    =  A   iln   2E^/A   -   v^  ^ ^-—^  -  A/w^.  (4.2) 

°  ^  12   A^  ^ 


To  find  spatially  localized  time  periodic  solutions  for  small  deviations 

?7  4 2 
from  the  dimerized  lattice  define-" '^^ 


A(x,t)  =  A^[l  +  5(x,t)]  (4.3) 


where 
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5(x,t)=e[A(X,T)exp  i-jt+A  (X  ,T)exp-l.jt]+£  5  (x,  t  )  +  e"5  (x,t ) .       (4.4) 


and  t  and  T  are  the  two  time  scales.  Using  multiple  time-scale 
asymptotic  perturbation  theory  one  finds  the  approximate  breather 
solution  (see  Appendix  C  for  details) 


5(x,t)  =  2e  sech(x/d'lcos  a)^t+-  e  seen'" (x/d)[ cos  2io„t-3l   (4.5) 

3  3  B   ■' 


where 


d  =  i;/2/2  e  (4.8a) 


and 


Wg  =  [l  -  ^  e"^)  Uj^.  (4.6b) 

?7  42 
It  13  easy  to  show-"'   that  the  effective  energy  density  is 

2A      '2     2        2    3    4 

^  -  ^^  1~T  ^  24  "^    '~^—   24'  ^^-^^ 

F    2^^ 

and,  thus,  that  the  classical  breather  energy  is 

E„  = -I  [^  -^^l^  -  oa\].  (4.8) 

B     IT  27 
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Breather  Dynamics  and  Optical  Absorption 

The  approximate  analytic  form  for  the  breather  (eq,  ^4.5)  was  used 
as  an  initial  condition  and  allowed  to  evolve  under  the  adiabatic  SSH 
equations  of  motion.   Figure  4.2(b)  shows  that  for  £=0.75  the  dynamics 
of  the  analytic  and  numerical  breathers  (fig.  4.2(a))  a'^e  very 
similar.   The  simpler  lattice  distortion  of  the  analytic  breather  may  be 
attributed  to  the  absence  of  higher  order  terms  in  e  in  the  analytic 
form,  eq.  4.5.   Equation  4.9  gives  the  classical  breather  energy 
as  -0.55  A  for  e=0.75  and  thus  essentially  accounts  for  the  energy 
"missing"  from  the  photoexcitation. 

Spatial  localization  of  nonlinear  excitations  implies  the  presence 
of  localized  electronic  levels  (in  one  dimension).  For  a  breather  there 
is  a  pair  of  intragap  levels  symmetrically  about  the  Fermi  level.  The 
lower  level  is  doubly  occupied,  the  upper  is  empty:  the  breather  is  a 
neutral  lattice  excitation.  From  fig.  4.3  it  is  clear  that  the  dynamics 
of  the  electronic  levels  for  the  numerical  and  analytic  breathers  are 
very  similar.   It  is  important  to  note  that  these  intragap  levels 
oscillate  deeply  into  the  gap  (for  this  numerical  breather  to 
-0.55  A  ).  However  the  intragap  levels  spend  much  of  the  breather  cycle 
close  to  or  beyond  the  band  edge,  thereby  increasing  the  effective  band 

gap. 

Transitions  into  and  out  of  these  localized  intragap  levels  produce 
contributions  to  the  optical  absorption  below  the  band  edge.   In  the 
absence  of  an  analytic  form  for  the  breather  electronic  wavefunctions  it 
is  not  possible  to  analytically  calculate  the  breather  optical 
absorption.   However,  useful  qualitative  insights  can  be  gained  from 
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analogies  with  the  exactly  calculated  (in  the  continuum  limit)  absorp- 
tion of  the  polaron.   The  breather  is  a  neutral  excitation,  whereas  the 
polaron  is  the  static  groundstate  of  the  system  with  a  single  charge 
added.   They  do,  however,  share  the  same  basic  topology  and  both  have 
localized  electronic  intragap  levels,  though  differing  in  their 
occupancies.   (The  relative  intensities  of  the  polaron  transitions  are 
shown  in  fig.  3.3-)  For   the  trans-polyacetylene  polaron  co  /A  -  0.7. 
For  the  breather  the  location  of  the  intragap  level  varies  with  time  and 
thus,  in  a  quasiclassical  sense,  an  average  in  cu  must  be  taken  over  a 

full  breather  period.   (Quantization  of  the  breather  levels  will  be 

(2) 
discussed  in  a  later  section.)  It  is  important  to  note  that  a        ,    the 

P 

analog  of  the  strongest  transition  in  the  polaron  is  absent  in  the 
breather  due  to  the  differing  electronic  occupancy.   For  this  numerical 

breather  the  intragap  level  moves  over  the  range  -0.55  -  1 .05  A  ,  over 

'  1  1 
which  a'  '  is  approximately  independent  of  the  location  of  the  intragap 

level.  The  relative  intensity  of  this  contribution  to  the  absorption 

is,  thus,  largely  governed  by  the  amount  of  time  spent  at  each  gap 

location.   Two  dominant  contributions  are  expected.   First,  a  strong 

contribution  close  to,  or  beyond,  the  band  edge,  where  the  intragap 

level  spends  most  of  its  period.   Second,  a  contribution  when  the 

intragap  level  is  farthest  from  the  band  edge  (for  this  numerical 

breather  at  -1.1  A  ).  There  should  be  a  weaker,  nearly  uniform 

(3) 
contribution  over  intermediate  energies.   The  transition  a    is  weak  at 

all  relevant  energies,  with  no  significant  contribution  at  the  lowest 

enersies.  For  this  breather  a  weak  contribution  at  -1.5  A  can  thus  be 
°  o 

expected.  The  nature  of  the  compensating  above  band  edge  bleaching  can 
also  be  understood  by  analogy  with  the  polaron.   If  the  polaron  is 
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narrow  it  bleaches  over  a  wide  range  of  energies — a  narrow  excitation 

requires  a  large  number  of  Fourier  modes,  including  those  at  high 

energies  and  wavevectors,  for  its  description.   4  wide  polaron,  on  the 

other  hand,  can  be  described  by  fewer,  lower  wave-vector  modes.  It 

bleaches  over  a  narrow  energy  range  close  to  the  band  edge.  The 

breather  is  extremely  wide  and  thus  can  be  described  by  only  the  low 

wave-vector  Fourier  components.   Its  above  band  edge  bleaching  should, 

therefore,  be  mainly  close  to  the  band  edge. 

Using  the  algorithm  described  in  Chapter  III  we  numerically 

calculated  the  optical  absorption  of  the  numerically  created  breather 

over  a  full  breather  period  (fig.  ^1.4).   Removal  of  the  breather  from 

the  environment  of  the  two  kinks  quantitatively  changes  the  kink 

dynamics  due  to  the  strong  modification  of  the  extended  wavef unctions  by 

the  kink.  The  expected  subband  absorption  at  ~A  can  be  seen.  There 

o 

are  also  contributions  close  to  the  band  edge,  which  overlap  with  the 
interband  absorption.   Here  the  sum  rule  is  obeyed  to  better  than  99.9i5 
when  averaged  over  a  full  breather  period,  and  to  ~99.6?  at  each 
instant.  Figure  4.5a  shows  the  change  in  optical  absorption  when  the 
numerical  breather  is  added  to  a  98  site  ring,  i.e.  fig  4.5a  =  fig.  4.4 
-  fig.  3.4.   The  breather  signature  is  clear:  absorption  enhancements  at 
-A  and  close  to  the  band  edge,  which  are  largely  compensated  by 
bleaching  of  the  interband  just  above  band  edge.    (It  should  be  noted 
that  the  optical  absorption  of  the  same  breather  in  the  presence  of 
a  KK  pair  is  a  single  peak  in  the  energy  range  1-2A  .  This  strong 
modification  is  due  to  the  effect  of  the  kinks  on  the  continuum  state 
with  which  the  breather  interacts.) 
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Figure  4.5b  shows  the  change  in  optical  absorption,  AaCul,  when  an 
analytic  breather  with  e=0.75  is  added  to  a  93  site  ring.   Below  the 
band  edge  is  a  single  broad  peak  arising  from  transitions  between  the 
lower  breather  level  and  the  conduction  band.   This  enhancement  is 
compensated  by  a  bleaching  over  a  narrow  energy  range  above  the  band 
edge. 

To  better  understand  the  dynamics  of  photoexcitation,  it  is 
interesting  to  follow  the  optical  absorption  on  the  subpicosecond 
timescale.  Figure  4.6  contrasts  the  time  average  of  the  optical 
absorption  over  the  periods  0.02  -0.04  psecs  and  0.30  -  0.34  psecs  after 
photoexcitation.  There  are  three  important  differences  to  note.   First, 
the  "midgap  levels"  are  moved  away  from  midgap  at  short  times — the  KK 
pair  is  still  evolving  and  the  intragap  levels  associated  with  them  have 
not  yet  fully  evolved.   Second,  because  of  this  splitting  there  are 
transitions  between  the  putative  midgap  levels.  This  is  characterized 
by  a  low  energy  contribution  to  the  absorption,  which  rapidly  decreases 
in  energy  and  intensity  as  the  intragap  levels  approach  midgap.  Third, 
and  perhaps  most  suprisingly,  the  intensity  of  transitions  between 
levels  which  will  eventually  evolve  to  the  breather  is  large  at  all 
times.  Thus,  the  electronic  properties  of  the  breather  fully  evolve 
very  rapidly  (<0.04  psecs),  although  the  characteristic  lattice 
distortion  of  the  breather  does  not  fully  evolve  until  t  -  0.15  psecs. 
(The  above  possibility  of  "coexisting"  breather  and  KK  complexes  is 
important  in  cis-like  materials  where  the  K  and  K  are  confined.) 
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Photoexci tation  in  the  Presence  of  Intrinsic  Gap  States 

To  this  point  only  photoexcitation  across  the  full  ►^and  gap  has 
been  considered.   However  ESR  and  electrical  conductivity-^-'^  experiments 
show  that,  even  for  "pristine"  samples  of  trans-polyacetylene,  there  is 
a  significant  density  of  spinless  charge  carrier3--kinks.   It  is 
therefore  interesting  to  study  photoexcitation  in  the  presence  of 
intrinsic  intragap  levels.   (Photoexcitation  in  the  presence  of 
extrinsic  intragap  levels  will  be  the  subject  of  Chapter  VI . 1 

First  consider  a  single  neutral  kink  on  a  chain.   Photoexcitation 
of  an  electron  from  the  midgap  level  to  the  bottom  of  the  conduction 
band  produces  the  electronic  configuration  of  a  negatively  charged 
polaron  and  a  positively  charged  kink.   The  rest  energy  of  the  polaron 
is  2/2  A  /tt  -  0.9  A  .   Photoexcitation  adds  energy  A   to  the  system. 
Thus  polaron  production  is  possible.   However,  for  a  pure  system  there 
is  symmetry  about  the  kink  and,  in  the  absence  of  an  alternative  seeding 
center  the  system  evolves  to  an  oscillatory  boundstate  of  the  kink  and 
polaron.   No  breather  is  produced. 

Now  consider  photoexcitation  of  an  electron  from  the  lower  to  the 

upper  intragap  level  of  the  polaron.   The  polaron  intragap  levels  are 

A  //2  from  midgap.   Thus  the  photoexcited  polaron  has  energy  -2.3  A 
0  o 

greater  than  the  groundstate.  Consequently  this  excited  polaron  can  be 
expected  to  be  unstable  to  the  formation  of  a  KK  pair  and  a  breather. 
Indeed  numerically  a  KK  pair  is  rapidly  produced  in  <0.05  psecs  and  a 
breather  is  also  created  (fig.  4.7).   By  energy  conservation  this 
breather  has  energy  -0.8  A  .  However,  a  simple  calculation  shows  that 
the  maximum  energy  of  the  classical  approximate  analytic  breather 
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is  -0.53  A  .   Thus  this  breather  is  not  well  described  by  eq.  '4.5, 

reflecting  the  need  for  higher  order  terms  in  e  expansion  for  the 

breather  (c.f.  results  for  the  breather  in,  e.g.  $   equation). 

Since  in  this  case  three  levels  oscillate  into  the  gap,  subband 

edge  absorption  can  be  expected  over  a  wide  energy  range.   Indeed,  fig. 

4.8  shows  that,  as  well  as  the  strong  midgap  absorption  from  the  two 

kinks,  there  is  a  low  energy  shoulder  from  kink-breather  transitions. 

There  is  also  an  enhancement  at  -1.5  A  arising  from  transitions  between 

o 

the  breather  and  continuum  levels.  The  intragap  bleaching  arise  from 
the  evolution  of  the  two  polaron  levels  to  midgap,  forming  the  KK  pair. 


Above  Band  Edge  Photoexcitation 

In  most  laboratory  experiments  photoexcitation  is  performed  with 
pumping  energy  well  above  band  edge,  thus  exciting  electrons  from  deep 
in  the  valence  band  to  high  in  the  conduction  band.  Such  highly  excited 
electrons  (and  holes)  can  be  expected  to  decay  towards  the  edges  of  the 
band  gap  by  optical  phonon  emission  with  a  time  constant  of,  typically, 
-  10    seconds.  However  close  to  the  band  edge  there  may  not  be  an 
optical  phonon  of  sufficiently  low  energy  to  mediate  further  decay. 
Further  decays  must  take  place  via  the  emisssion  of  acoustic  phonons , 

_Q 

which  typically  have  a  time  constant  of  -10    sees.   In  the  absence  of 
a  preferred  non-radiative  decay  route,  via  e.g  coherent  electron-hole 
scattering,  this  "bottleneck"  may  mean  that  the  system  stays  in  a 
metastable  excited  state  that  can  then  evolve  adiabatically  to  form 
nonlinear  excitations.   Amongst  these  nonlinear  excitations  one  may 
expect  to  find  KK  pairs,  polarons,  polaron-polaron  bound  states. 
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breathers  and  excltons.   The  choice  of  channel  will  depend  on,  e.g. 
correlation  and  impurity  effects. 

A  f ally  quantum  non-adiabatic  model  is  not  yet  available  to 
investigate  this  problem  in  detail.   Yet  useful  insights  can  be  gained 
from  a  simple  model,  which  mimics  some  of  the  nonadiabatic  effects.  At 
t=0  an  electron  was  excited  vertically  (Ak=0)  from  the  fourth  highest 
level  in  the  valence  band  to  the  fourth  lowest  level  in  the  conduction 
band.   The  excited  electron  and  hole  were  then  allowed  to  decay  towards 
the  band  edge  independently  and  randomly  with  time  constant  t. 
Baeriswyl  -"  has  estimated  that  for  trans-polyacetylene  t  ~  0.015 
psecs.  We  simulated  only  the  simplest  model  in  which  the  electron  and 
hole  can  decay  rapidly  all  the  way  to  the  band  edge.  Figure  4.9  shows 
that  when  the  electron  and  hole  are  excited  the  system  tries  to 
equilibriate  by  changing  its  dimerization.   In  this  simulation  the 
electron  and  hole  had  both  reached  the  band  edge  by  0.1  psecs  after 
which  a  KK  pair  were  rapidly  produced. 

Neutron  Scattering  Cross-Section  of  the  Breather 

Inelastic  scattering  of  low  energy  neutrons  should  provide  a  direct 
probe  of  the  lattice  structure  of  polyacetylene. 

In  this  section  a  simple  phenomenology  is  used  to  elucidate  the 
basic  features  of  the  classical  dynamic  structure  factor  of  the 
breather.  It  is  assumed  that  any  lattice  distortion  can  be  considered 
as  a  superposition  of  independent  collections  of  phonons,  kinks, 
polarons,  breathers  etc.  Although  the  focus  here  is  on  breathers, 
similar  phenomenologies  can  be  developed  for  the  kink  and  the  polaron. 
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In  both  cases  there  are  contributions  to  the  struct jre  factor  at  oj  ~  0, 
as  well  as  responses  associated  with  the  internal  vibratory  modes. 
Here,  as  has  been  done  previously  for  the  sine-Gordon  system,   we  treat 
the  collection  of  breathers  as  an  ideal  gas,  in  which  all  the  internal 
structure  is  reflected  in  a  q-dependent  form  factor.   An  outline  of  the 
derivation  is  given  here;  full  details  are  in  Appendix  D. 
Define  the  classical  structure  factor  as 


S(q,a))  =  —   J  dx  dx  dt  dt  exp[ia)(t  -t  )-iq(x  -X  )] 

(2tt) 

(4.9) 

<  exp  iqu(x  t  ■)  exp-iqu(x^t  )> 


where  u  is  the  lattice  displacement  at  location  x  and  time  t.  Keeping 
only  one  phonon  processes  and  assuming  that  the  thermal  average  can  be 
implemented  by  integrating  over  a  classical  distribution  P(v^,  then 


_-,,,  n    (o)    ;T)  _ 

S(q,(jj)    =   e   " r —   [dvP(v)q   [jdx.dt,    exp   iY    (q+Tr/a-v/jo/c    ") 

(2Tr ) 


(4.10) 


exp   iY    {u>-qv)t^[u{x^t^)]][x^^x^,t^->-t^] 


-2W 
where  e    is  the  Debye-Waller  factor,  c^  is  the  velocity  scale  and 

u(x,t)=na+(-1 )"  u(x,t).   For  simplicity  we  use  the  Maxwell-Boltzmann 

distribution,  p'  (v).   After  a  little  algebra  one  finds 
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(27r)^  I  q  1  q 


(i^.n) 


*fJiif'(^;)iv^^[^)H-^(^-)iv^^[:^).36|f^(.,)ip"^./q)i] 


wnere 


-1  (jjtri'j) 

k^  =   d   Mq   -   -/a   -^  (-— ^)]                          (4.12a) 

c  ^ 
o 


k     =   d   ^    [q  +   ir/a   -  w^/c^^q]  (4.12b) 


and 


f^(q)   =   TTd   sech    (iTqd/2)  (4.13a) 


f    (q)    =   irqd^   cosech    (TTqd/2)  (4.13b) 


For  simplicity  we  look  at  the  behavior  of  S(q,(ij)  in  the  T=0  limit, 
i.e.  when  P^   (v)  ->■  5(v).   Taking  this  limit  is  not,  of  course,  strictly- 
meaningful  for  a  classical  calculation  but  does  help  display  the 
essential  features  of  equation  4.11. 

We  find 


S(q,uj=0)/S(q,a)=±2a)„)  =  3&.  (4.14) 

o 


This  zero  frequency  inelastic  scattering  can  be  observed  separately  from 
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the  elastic  scattering  of  the  phonons .   Further  contributions  at  co=0  are 
expected  from  other  nonlinear  excitations.   These,  however,  may  have 
differing  temperature  and  intensity  dependences. 
Also 


S(q,u)=±u)„)/S(q,io=±2a)„)  =  ^      /^  tanh^[''^^'^-(  1 -q/ir )  ]        (4.15) 


For  nearly  all  values  of  q  (except  q  ~  it)  this  gives  almost 
independently  of  e 


S[q,iij  =  ±w„]/s(q  =  ±T,  to  =  ±2a)„]  ~  9  x  10  ^       {^.^6) 

D  D 


At   the   zone   boundary   q=±iT   one  finds 


S[q   =   ±T,   10  =   ±u„)/S(q   =   ±t  ,   y  =   ±2u)„]    =  -^  -    3-8        (4.17) 

25/e 


Quantization  of  the  Breather 

One  of  the  major  assumptions  of  the  work  has  been  that  the  lattice 
degrees  of  freedom  can  be  treated  as  classical  variables.   One  of  the 
effects  of  a  quantum  treatment  is  to  renormalize  the  mass  of  the 
nonlinear  excitations:  e.g.  the  mass  of  the  kink  is  reduced  by  ~  25?  by 
quantum  effects.^   In  addition  to  the  lattice  being  quantized  the 
locations  of  the  electronic  levels  are  quantized,  lying  within  the 
bounds  of  the  range  of  values  of  their  classical  analogs.  Due  to  the 
electron  hole  symmetry  of  the  formalism,  the  kink  intragap  level  is  at 
midgap  for  both  the  classical  and  quantum  treatments.   The  locations  of 
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the  classical  extended  states  and  the  intragap  polaron  levels  are 
essentially  time  independent  and  thus  the  locations  of  the  quantized 
levels  are  well  defined.   For  the  breather  the  situation  is  entirely- 
different:   it  is  an  inherently  dynamic  excitation,  whose  lattice 
distortion  and  electronic  spectrum  are  time  dependent. 

In  a  previous  section  we  derived  an  approximate  analytic  form  for 
the  energy  of  the  breather.  This  can  easily  be  quantized  at  the  Bohr- 
Somraerfeld  level  by  demanding  that 


J  =  J  p  dA  =  2mT  (4.18) 


where  J  is  the  action  and  p   is  the  conjugate  momentum  to  the  order 
parameter,  A.  Thus 


2  r  A^ 

2mr  =  — -  J  ^  dx  dt  (4.19) 

TTV  ^     d 

F   ^^ 


which  gives 


2nTr  =  —  i  [l  -  ^  l^\  (4.20) 


-R        27 


Inverting  and  substituting  into  eq.  4.8  for  the  breather  energy  gives 


2  2  2 

n  TT  O)  . 

E  =  nto   1  -  ^  +  0(0")]  (4.21) 

^       72A  2  , 
o 

Equation  (4.21)  shows  explicitly  that  a  breather  of  quantum  number  n  is 
energetically  favored  over  n  incoherent  phonons. 
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To  calculate  the  optical  absorption  spectrum  of  the  quantized 
breather  it  is  necessary  to  calculate  the  energy  of  the  excited 

U   1 

breather.   One  might  hope  to  do  this  using  an  effective  Lagrangian   of 

the  type  used  for  the  groundstate  breather:   there  are  however  technical 

difficulties  in  extending  this  to  excited  states.   Other  approaches 

might  be  a  generating  functional  formalism  ■''  °  or  by  direct  numerical 

quantization  of  equation  4.19.  Even  if  this  can  be  achieved  there 

remains  the  conceptual  problem  of  quantizing  the  electronic  spectrum, 

from  a  knowledge  of  the  breather  energy  spectrum. 

At  present  we  can  only  say  that  on  quantization,  the  location  of 

the  breather  electronic  levels  and  thence  the  energy  of  contributions  to 

the  optical  absorption  spectrum  are  bounded  by  the  limits  of  their 

classical  analogues.  Thus  for  the  quantized  breather  we  expect  a  single 

sharp  absorption  peak  somewhere  in  the  range  -  1-2  A  . 

o 

Comparison  between  Theory  and  Experiment 

In  this  chapter  a  simple  scenario  has  been  proposed  for  the  short 
time  evolution  of  the  photoexcited  system.  Before  proceeding  to  further 
detailed  studies  of  the  dynamics  in  trans-polyacetylene  and  related 
materials  it  is  important  to  ascertain  what,  if  any,  experimental 
evidence  there  is  to  support  the  validity  of  this  scenario. 

Figure  i| .  1 0  shows  the  experimentally  measured  change  in  absorption 
during  photoexcitation  of  trans-polyacetylene  at  1 0K .  ' '   The  spectrum 
shows  three  basic  features:  a  broad  peak  at  ~0.45eV,  a  narrow  peak 
at  -1.35eV  and  bleaching  over  a  wide  energy  range  above  the  band  edge. 
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Clearly  we  rnjst  identify  the  O.^SeV  line  with  the  kink  "rnidgap" 
absorption  and  the  1.35eV  line  with  the  breather. 

The  O.Jj5eV  Line.^^ 

As  well  as  the  above  mentioned  high  energy  peaks,  a  large  n^jmber  of 
i.r.  active  modes  have  been  observed  during  photoexcitation.  i»^o»3^  j^^ 
experiments,  in  which  the  temperature,  laser  intensity  and  chopping 
frequency^  have  been  varied  some  of  these  peaks  have  been  shown  to  be 
directly  correlated  with  the  O.^SeV  peak.  This  demonstrates  that  they 
arise  from  the  same  center,  and  that  the  center  is  charged.  Further 
certain  of  these  i.r.  modes  have  been  convincingly  identified,  in 
theoretical  studies  by  Horovitz,'^"  as  arising  from  the  pinned  modes  of 
an  excitation,  which  breaks  the  lattice  translational  symmetry.  This 
does  not  uniquely  identify  the  excitation  as  a  kink:  it  might  equally 
well  arise  from  a  polaron.   Recently  Horovitz's  work  has  been  extended^^ 
and  it  has  been  shown  that  there  are  additional  contributions  to  Lhe 
i.r.  absorption  unique  to  the  kink,  arising  from  its  internal  modes. 
This  may  allow  for  the  unique  identification  of  the  origin  of  the  0.45eV 
line  from  optical  data  alone.  From  the  intensity  ratios  of  certain  of 
the  i.r.  lines,  the  mass  of  the  charge  center  has  been  estimated 
at  -lOmg.  This  is  in  good  agreement  with  estimates  for  both  the  kink 
and  polaron  masses. 

The  determination  of  the  spin  of  the  charge  carriers  will 
differentiate  between  the  polaron  and  the  kink.   Infrared  studies  have 
seen  the  same  features  in  doping  as  in  photoexcitation. 35  ■p^e  features 
are  independent  of  the  kind  of  dopant  and  thus  must  arise  from  intrinsic 
excitations  of  the  system.  Further,  on  doping  a  midgap  absorption  peak 
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appears  correlated  with  some  of  the  i.r.  peaks  associated  with  the 
midgap  absorption  produced  during  photoexcitation.   This  shows  that  the 
spectroscopic  features  arising  during  photoexcitation  and  on  doping  are 
associated  with  the  same  lattice  excitation.   Flood  and  Heeger  showed 
in  ESR  experiments  that  the  number  of  spin  carriers  does  not  increase 
during  doping.   This  strongly  supports  the  identification  of  the  0.^5  eV 
line  as  arising  from  charged  spinless  kinks. 

The  1  .35eV  Line.'^^ 

No  i.r,  active  mode  has  been  found  to  be  correlated  with  the  1.35eV 
line,  implying  that  this  peak  arises  from  a  neutral  excitation.   Sub- 
picosecond  spectroscopy  shows  that  the  1.35eV  line  is  strong  and  already 
decaying  0.1  psecs  after  photoexcitation.-'  -'^   This  short  timescale  for 
production  and  the  estimate  of  the  quantum  production  efficiency  as 
being  of  order  unity  strongly  suggest  that  it  arises  from  an  intrinsic 
excitation  of  the  lattice.  All  of  the  above  are  in  qualitative  agree- 
ment with  the  breather  model:  the  breather  is  indeed  neutral  and  its 
electronic  signature  is  well  developed  after  one  phonon  period  (~0.04 
psecs).  The  observed  intensity  decrease  may  be  attributable  to  energy 
loss  to  e.g.  the  vibronic  modes  of  the  C-H  bond.   There  are  four  other 
important  pieces  of  data.   First,  the  1.35eV  line  decreases  in  intensity 
with  increasing  temperature  and  entirely  disappears  at  ~  150K.  ' 
Second,  at  approximately  the  same  temperature  there  is  a  sharp  increase 
in  the  photocurrent  .-'^   Third,  the  energy  of  the  peak  decreases  by 
approximately  3%  on.   deuteration.^"  Fourth,  the  excitation  responsible 

ilQ 

for  the  1.35eV  line  is  dipole  forbidden  from  the  groundstate.  '^ 


76 


Let  us  examine  each  of  these  in  turn.   First,  the  breather  is  a 

lattice  excitation  and  thus  one  expects  it  to  be  strongly  temperature 

dependent.   Indeed  the  temperature  at  which  the  breather  becomes 

unstable  can  be  estimated  by  a  simple  argument.  A  typical  breather  may 

have  energy  -0.5  A   (=0.35eV)  and  be  localized  over  -2^=15  lattice 
o 

sites.  Making  the  plausible  assumption  that  the  breather  becomes 

unstable  when  the  thermal  energy  is  equal  to  the  breather  energy  gives 

an  estimate  of  -3OOK  for  the  temperature  at  which  the  breather  becomes 

unstable.   Second,  the  increase  in  the  photocurrent  at  -I5OK  cannot  be 

understood  within  the  breather  model.  Although  it  has  not  been  shown 

that  the  onset  of  the  photocurrent  and  the  disappearence  of  the  breather 

are  correlated,  this  must  be  regarded  as  a  failing  of  the  breather 

model.   Third,  assuming  that  the  breather  in  the  hydrogenated  and 

deuterated  systems  has  the  same  quantum  number,  one  expects  the  energy 

of  the  peak  to  decrease  by  [  1-M(CH)/M(CD)  ]  ' '^^  -  1^%.      Fourth,  the 

breather  has  the  electronic  configuration  of  the  groundstate  and  can 

decay  to  it  only  by  phonon  emission. 

A  quantitative  comparison  between  the  experimentally  observed  and 

numerically  calculated  optical  absorption  data  is  difficult  for  a  number 

of  reasons:  (a)  the  presence  of  electron-electron  interactions  shifts 

the  "midgap"  absorption  to  -0.45eV;  (b)  in  the  absence  of  a  full  quantum 

theory  of  the  breather  one  can  only  say  that  the  breather  absorption 

should  consist  of  a  sharp  peak  at  an  energy  within  the  bonds  of  its 

classical  energy;  i.e.,  1.0-2.10  A  ;  (c)  we  have  assumed  a  breather 

o 

density  of  ~]%:    the  actual  breather  density  is  unknown;  (d)  detailed 
structure  arising  from  vibronic  modes  of  the  C-H  unit  are  not  included 
in  this  model. 
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Although  a  quantitative  comparison  is  not  possible  the  qualitative 
consistency  of  the  kink-breather  photoexcitation  scenario  with  data  of 
many  different  types  clearly  justifies  further  study  of  this  model.   In 
the  next  chapter  we  discuss  the  modifications  to  the  S3H  Hamiltonian 
needed  to  study  cis-polyacetylene .  We  then  explore  the  dynamics  of 
photoexcitation  in  cis-polyacetylene. 
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Figure   4.4 

Average  optical  absorption,  ciCu),  over  one  phonon  period  and 
over  a  98  length  chain  for  the  numerical  breather.  For  the 
dimerized  lattice  see  fig  3.4. 
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Figure   4,5 

Change  in  optical  absorption  when  (a)  the  numerical  breather 
(solid  line)  or  (b)  the  analytic  breather  (dashed  line)  is  added  to 

the  98  site  ring. 
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Figure  H.6 

Change  in  optical  absorption,  AaCw),  on  photoexcitation  for 
trans-pol'yacetylene  for  (a)  t=21-43  fsecs  (solid  line) 
(b)  t=301-343  fsecs  (dashed  line) 
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F  i  gur  e  4.8 

The  change,  Aa(w),  in  the  optical  absorption  on 
photoexcitation  of  the  electron  polaron. 


Figure  4.9 

Above  Band  Edge  Excitation 

(a)  Dynamics  of  system  with  electron  hole  pair  excited  above 

band  edge. 

(b)  Location,  with  respect  to  the  Fermi  level,  of  electron  and 

hole  (solid  line).  Level  occupied  by  photoexclted  electron 

and  hole  (dashed  line),  where  the  band  edge  state  are  levels  49  and  50, 
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Figure  4.10 

Experimentally  observed  change  in  absorption  on  photoexcitation  of 
(a)  Gis-polyactylene  and  (b)  trans-polyacetylene  at  10  K. 
From  data  of  Orenstein  et  al .  ''^^ 


CHAPTER  V 
DYNAMICS  IN  CIS-POLYACETYLENE  AND  RELATED  MATERIALS 


Unlike  trans-polyacetylene  the  two  isomers  of  cis-polyacetylene  are 
not  degenerate;  the  trans-cisoid  configuration  being  energetically 
preferred  over  the  cis-transoid  configuration  (fig.  5.1).  This  non- 
degeneracy  can  be  viewed  as  arising  from  the  differing  orientation  of 
the  C-H  bonds  with  respect  to  the  C-C  bonds  in  the  two  isomers. 
Alternatively  one  may  view  the  non-degeneracy  as  arising  from  the 
differing  third-nearest  neighbor  interactions.   As  the  trans-cisoid 
configuration  has  6E  energy  per  site  less  than  the  cis-transoid 
configuration  a  KK  pair  at  a  distance  i   apart  are  confined  by  an 
energy  5E.5,.  Single  kinks  may,  however,  still  stabilize  at  chain 
ends.   Polarons  and  breathers,  which  do  not  involve  a  change  of  phase 
are  still  stable. 


Statics  in  Cis-Polyacetylene 


The  non-degeneracy  of  cis-transoid  and  trans-cisoid  can  be  modeled, 
in  the  continuum  limit  by  the  addition  of  an  extrinsic  gap,*^  A  ,  to  the 
intrinsic  gap  produced  by  the  electron-phonon  coupling.  This  gives  the 
total  gap,  A(x)  ,  as 


A(x)  =  A^(x)  +  Ag.  (5.1) 
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This  does  not  alter  the  structjre  of  the  electronic  spectram  or  it: 
eigenf unctions .   It  does,  however,  modify  the  self-consistency  gap 
relation  to 


2 

A.(x)  =  A(x)  -  A   =  -  -^   E'  [j*(x)  v(x)  +  v\x)  u(x)].  (5.2) 
1  e       2      ^  ■' 

wQ   v,s 


Substitution  of  the  eigenf unctions  of  the  purely  diraerized  lattice  (eq. 
2.17)  gives  the  implicit  relation  for  A  , 


A^  =  XA^  m  {^/^).  (5.3) 

e     o      00 
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It  can  easily  be  shown^  rigorously  (or  by  direct  substitution)  that  the 

kink  is  no  longer  a  solution  of  eq.  5.2.   One  now  finds  that,  as  well  as 

the  polaron  solution  to  eq.  5.2  there  are  also  bipolaron  solutions,  with 

either  zero  or  two  electrons  in  each  intragap  state.   There  is  also  an 

"exciton"  solution  with  a  single  electron  in  each  intragap  level.   In 

trans-polyacetylene  both  the  bipolaron  and  the  exciton  are  unstable  to  a 

kink-antikink  pair. 

In  the  SSH  model  the  necessary  symmetry  breaking  can  be  modeled  by 

a  term,  which  in  the  continuum  limit,  reduces  to  the  extrinsic  gap. 

Here  two  possible  additions  to  the  SSH  Hamiltonian  are  considered 


"l  =  S  ^  (-^^"  tCl  'n'<^-n.^^  ^5.^) 


and 


"3  =  S   ^   K.^^'n^^^n^,^'  ^5.5) 


n.3- 


even 
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It  is  straightforward  to  show  that  in  the  continuum  limit 


t^  =  2Ag  =  2t2.  (5.6) 


Although  both  Hi  and  H-,  modify  the  continuum  equations  in  an  identical 
manner,  they  may  behave  qualitatively  differently  in  the  discrete 
system.  In  fig.  5.2  the  potential  energy  (electronic  plus  lattice)  is 
plotted  as  a  function  of  the  dimerization  for  H^  and  Ho.  The  parameters 
t-i  and  to  were  chosen  to  numerically  give  the  same  total  groundstate 
gap.   Close  to  the  absolute  minimum  the  two  potentials  are  very  similar 
and  the  dynamics  within  them  should  be  the  same.  Although  the  potential 
around  the  metastable  minimum  is  different  in  the  two  cases,  this  is 
less  important  for  the  dynamics  of  confined  "kinks".   As  the  relevant 
dynamics  should  be  qualitatively  similar  for  the  two  models  only  Hi  will 
be  considered  here. 

It  has  not  yet  been  possible  to  derive  unique  values  for  the 
parameters  of  cis-polyacetylene  from  the  experimental  data.   In  this 
study,  therefore,  the  rescaled  parameters  of  trans-polyacetylene  are 
used.   The  dynamics  are  then  studied  as  a  function  of  the  confinement 
strength  r  =  t  /t  .   (From  resonant  Raman  scattering  studies  Vardeny  et 
al .   have  estimated,  assuming  that  the  electron-phonon  coupling 
constant  (X)  is  the  same  for  cis-polyacetylene  as  for  trans- 
polyacetylene,  that  the  extrinsic  gap  is  -  5?  of  the  total  gap.  Such  a 
gap  increase  corresponds  to  r  -  0.01,  for  the  parameter  set  used  here. 
As  we  shall  see  this  suggests  that  the  kinks  in  cis-polyacetylene  may  be 
weakly  confined.) 
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Dynamics  in  Cis-Polyacetylene 

The  potential  energy  of  the  dimerized  lattice  was  found  to  vary 
linearly  with  r  over  the  range  of  values  studied  (0<  r  <  CDS'). 
(Although  this  appears  to  be  a  narrow  range  of  values,  we  shall  see  that 
for  r=0.06  the  kinks  are  strongly  confined.)  As  free  kink  production  is 
forbidden  we  expect  that  photoexcitation  will  lead  to  the  production  of 
a  KK  bound  state.  We  found  that  there  are  two  regimes  in  r,  with 
qualitatively  different  behavior  during  photoexcitation: 

1.  In  the  "weak  confinement  regime"  (r  <  0.015)  a  KK  pair  is 
formed.  The  kinks  separate  (fig.  5.3a)  to  some  maximum 
separation.   They  then  evolve  to  a  dynamic  cycle  of  quasi- 
independent  kink  and  antikink  which,  after  a  long  time  is 
expected  to  produce  either  (a)  a  dynamic  complex  involving 
kink-antikink  creation  and  recombination,  or  (b)  a  static 
configuration  with  the  kink  and  antikink  at  some  finite 
equilibrium  separation. 

2.  In  the  "strong  confinement  regime"  (r  >  0.015)  it  is  not 
possible  to  separately  identify  a  kink  and  an  antikink  (fig. 
5.3b):  the  system  evolves  to  a  single  complex  lattice 
excitation.  The  electronic  signature  of  this  complex  is  a  pair 
of  singly  occupied  intragap  states,  that  oscillate  over  a 
narrow  energy  range  around  ~A   from  midgap.   This  is  clearly 
the  exciton,  predicted  by  the  continuum  theory.  What  is  more 
suprising  however,  is  that  the  new  band  edge  states  also 
oscillate  periodically  into  the  gap.   This  breather-like 
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electronic  signature  cannot  be  identified  with  a  specific 
lattice  excitation,  bit  appears  to  be  an  integral  part  of  the 
photoexcited  complex. 
The  optical  absorption  was  calculated  for  the  above  two  examples. 
Figure  5.4  shows  the  difference,  Aado),  between  the  long  time  evolved 
photoexcited  system  and  the  groundstate.   For  r=0.002  the  absorption 
(fig.  5.4a)  is  qualitatively  similar  to  that  of  trans-polyacetylene 
(fig.  4.6bl.   The  major   difference  is  that  the  kink  midgap  and 
therefore,  compensating  above  band  edge  absorption  is  strongly 
diminished.   This  makes  the  high  energy  shoulder,  arising  from  the 
breather,  much  more  noticeable.   It  is  clear  that  for  r=0.06  the 
absorption  spectrum  is  qualitatively  different.   First  the  high  value  of 
r  increases  the  band  gap  to  -2.5  ^    .      Second  the  contributions 
at  -  A  /2  and  -  5A  /4  arise  from  transitions  between  the  putative  midgap 
levels  and  between  the  midgap  levels  and  the  continuum  states 
respectively.   The  broad  peak  at  -1.5-2.0  A   arises  from  transitions 
within  the  breather-like  modes  of  the  photoexcited  system. 

We  have  seen  numerically  that  in  cis-polyacetylene  there  is  an 
excitation  analogous  to  the  breather  of  trans-polyactylene.   Here  for 
high  r,  however,  it  is  part  of  an  "exciton-breather"  complex.   As  the 
breather  is  coupled  to  an  electronic  excitation  we  do  not  expect  its 
optical  signature  to  be  strongly  temperature  dependent.   Indeed,  in  the 
real  system  (fig.  4.10a)  there  is  a  large,  temperature  independent 
peak  '  '   just  below  the  band  edge,  which  has  previously,  with  little 
experimental  support,  been  ascribed  to  an  excitonic  breather.   In  this 
simulation  we  also  see  contributions  to  the  optical  absorption  close  to 
midgap.  There  is  no  analog  of  these  in  the  experimentally  observed 
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optical  spectra.   (There  are  small  enhancements  in  the  absorption  close 
to  midgap  bJt  they  are  due  to  remnants  of  trans-polyacetylene  in  the 
system.')  This  qualitative  disagreement  may  be  due  to  a  bad  choice  of 
parameters  for  cis-polyacetylene:   with  some  other  suitable  choice  the 
dynamics  may  be  different.   It  is  also  possible  that  the  Hamiltonian 
used  to  describe  the  symmetry  breaking  is  inappropriate.   It  is, 
however,  possible  that  the  scenario  described  above  is  qualitatively 
incorrect  for  this  system. 

Breather  Dynamics  in  Cis-Polyacetylene 

In  fig.  5.5  the  period  of  an  £=0.75  breather  is  plotted  as  a 
function  of  r,  for  breathers  in  both  the  absolute  and  metastable 
minima.  As  r  increases  the  potential  around  the  absolute  minimum 
becomes  deeper  and  narrower,  thereby  decreasing  the  period  of  the 
breather.   The  metastable  minimum  becomes  shallower  and  wider  with 
increasing  r.  The  period  of  the  breather,  therefore,  increases.   Indeed 
for  r>0.05  the  metastable  minimum  becomes  so  shallow  that  it  is  no 
longer  able  to  support  the  breather. 

Dynamics  in  Finite  Polyenes 


The  finite  polyenes   ((CH)^^,  n=2,H  ,6 , . . .)   are  the  short  chain 
analogs  of  polyacetylene.  As  the  polyenes  (and  polyacetylene)  belong  t 
the  Cpf,  symmetry  group  their  single  particle  spectrum  consists  of 
alternating  'a_  and  ^B.,  states.  Since  in  the  groundstate  all  of  the 
orbitals  are  doubly  occupied  the  groundstate  should  have  'A  symmetry 
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and  the  first  excited  state  should  have  3,^  symmetry  with  a  dipole 

allowed  transition  between  them.  In  studies  of  octatraene  it  has  been 

shown  that  the  optical  absorption  band  edge  is  above  the  band  edge  for 

fluorescence.   These,  and  other  data,  strongly  suggest  that  there  is  a 

reversal  of  level  ordering  and  that  the  first  excited  state  is  A^,  and 

g 

thus  dipole  forbidden  from  the  groundstate.   In  photoexci tation  the  next 
highest  level,  which  has  B.^  symmetry  is  excited.   This  promptly  decays 
by  a  radiationless  route  to  the  excited  A^.  state  which  then  fluoresces 
to  the  groundstate.  This  reversal  of  levels  cannot  be  explained  in  the 
single  electron  theory,  but  is  well  explained  by  theories  which  include 
electron  correlations  e.g.  the  Pariser-Parr-Pople  theory. 

In  this  study  polyacetylene  has  been  viewed  exclusively  in  the 
single  particle  picture.   As  mentioned  in  Chapter  II  this  approximation 
is  satisfactory  for  the  understanding  of,  at  least  some  of,  the 
experimental  data.  Of  course  it  would  be  more  satisfactory  if  the 
polyenes  and  polyacetylene  could  be  understood  within  the  same 
theoretical  framework.   Although  the  natural,  and  qualitatively 
successful,  approach  to  the  polyenes  remains  the  strongly  correlated 
models,  it  is  instructive  to  consider  the  nonlinear  dynamics  of  the 
single  particle  theory  for  the  polyenes.   Although  it  is  unlikely  that 
nonlinear  dynamics  alone  can  adequately  describe  polyene,  this  study  may 
elucidate  the  role  nonlinear  excitations  in  the  short  chain  system. 
There  are  two  competing  effects:  on  the  one  hand  a  short  system  has 
fewer  phonon  modes  to  act  as  energy  sinks  for  energy  removed  from 
nonlinear  excitations;  on  the  other  hand  in  a  system  with  so  few  modes 
the  idea  of  a  nonlinear  lattice  excitation  may  not  be  meaningful. 
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The  Dlmerized  Lattice 

Figure  5.6  shows  that  the  lattice  energy  per  bond,  E(N),  varies 
inversely  with  the  chain  length;  and  that  the  configuration  with  a 
double  bond  at  each  chain  end  is  lower  in  energy.  Thus  the  finite 
polyene  is  intrinsically  cis-like.   Indeed  for  N<20  the  configuration 
with  single  bonds  at  each  chain  end  is  absolutely  unstable  to  the 
configuration  with  double  bonds  at  each  end.   As  in  cis-polyacetylene , 
therefore,  single  kinks  are  not  free  and  single  kink  dynamics  are  not 
expected  to  be  important. 

Breather  Dynamics  in  Finite  Polyenes 

The  lifting  of  the  groundstate  degeneracy  does  not  qualitatively 
affect  the  breather.  A  numerical  study  of  the  dynamics  of  a  £  =0.75 
analytic  breather  on  a  chain  of  10  atoms  (with  SSH  parameters)  showed 
that  the  breather  is  persistent  to  at  least  1  psec.  Figure  5.7  shows 
the  difference  between  the  optical  absorption  of  the  chain  with  and 
without  a  breather.  It  is  evident  that  the  optical  absorption  is  very 
similar  to  that  of  trans-polyacetylene — a  single  subband  edge  absorption 
peak  compensated  by  an  above  band  edge  bleaching. 

Photoexcitation  in  the  Presence  of  Damping 

Although  emission  of  optic  and  acoustic  phonons  from  a  breather  during 
photoexcitation  is  small,  non-radiative  decay  routes  via  e.g.  vibronic 
modes  of  the  C-H  bond  or  a  carbon-impurity  bond  may  provide  an  efficient 
sink  of  energy.   Unfortunately,  due  to  difficulties  in  the 
interpretation  of  the  experimentally  measured  phonon  dispersion 
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curves, '^-^  no  estimate  of  the  strength  of  dissipation  in  polyacetylene 
has  been  made.   As  a  crude  approximation  we  model  all  forms  of 
dissipation  by  the  addition  of  a  spatially  uniform  velocity  dependent 
damping  term  to  the  SSH  equations  of  motion  (eq  3.11).   This  is,  of 
course,  highly  phenomenological  as  dissipative  effects  are  probably  only 
spatially  uniform  over  macroscopic  length  scales.   Further,  rather  than 
simply  providing  an  energy  sink,  vibronic  modes  of  the  C-H  bond  may  add 
a  crystal  phonon  field  to  the  problem. 

Figure  5.8  shows  the  result  of  single  electron  photoexcitation  in 
the  presence  of  damping  of  strength  u  =  0.15  w  .   A  kink  and  antikink 
are  formed.   Breather  production  is  suppressed. 
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Schematic  of  lattice  -  upper  is  trans-cisoid,  lower  is  cis-transoid , 

Total  potential  energy  of  cis-polyacetylene.  Note  that  the 

trans-cisoid  isomer  is  lower  in  energy. 


CD 


cu 

x: 

1 

„ 

4-> 

CO 

c 

CO 
CU 

c 

0) 
X3 

o 

L 

-o 

c 

cu 

CO 

c 

c 

o 

o 

0 

• 

<4_, 

'J) 

cu 

c 

c 

CO 

0) 

c 

•o 

J-> 

1 

J-0 

ra 

_ 

cu 
cu 

■3 

cu 
o 

•  i-H 

or 

CO 

1. 

>, 

V 

.— 1 

s 

cu 

CO 

c 

•^ 

c 

ex 

T3 

0) 

to 

1 

r— f 

CO 

cu 

>, 

-.— ( 

o 

cu 

■o 

o 

4-3 

o 

CO 

Cm 

J-) 

CO 

c 

CX3 

>. 

,— 1 

c 

c 

0 

o 

CU 

00 

c 

O 

Q. 
1 

c 

CO 

4-> 

c 

•n4 

to 

o 

4-> 

o 

X3 

r— 1 

o 

C 

1 — 1 
0^ 

C 

c 

CO 

Q 

3 

C^l 

03 

CU 

c 

1— f 
CO 

o 

C 

3 

•f-J 

a-3 

CO 

.—1 

4-> 

CO 

T3 

2 

cC 

H 

-P 

•  r-( 

O 

C. 

JJ 

CO 

' 

o 

c 

— 

CtH 

CU 

CU 

c 

cu 

C 

c 

c 

Q. 

cu 

,-\ 

CO 

O 

CU 
Q. 

-U 

cu 
o 

CO 

cu 

4-> 

o 
■o 

CO 

>> 

rH 

OO 

C 

1 

(>-i 

X 

105 


\ 


.\ 


/ 


/ 


m 

V 

rs 


C 
0) 

r-{ 
>, 
■i-> 

<D 
O 

CC 

>, 

r-H 
C      OJ 

a  o 
I  o 
fO      • 

•rH     o 

o 
II 

o  t- 

c  — 

O     CO 

4-> 

CD 

•rH 

o 

X 

0) 

o 
o 


107 


AX 


3. 


CD 

5 

bO 


0) 

c 

I— I 
>> 

o 
ta 

0  • 
Q.  vo 

1  O 

m  • 
■^  o 
o 

II 
J- 
o  t- 

c  ^■ 

O    J3 

•r^  ^--' 
J-> 

CC 
JJ 

•iH 
O 

X 

cu 
o 

o 

CL, 


109 


o 


AX 


3^ 


no 


40   - 


3 


0   = 


-40   - 


-80   - 


2.0 


4.0 


6.0 


E/A, 


Figure  5.^ 

Change  in  optical  absorption,  AaCwl,  on  photoexcitation  for 
r  =  0.002  (solid  line)  and  r  =  0.06  (dashed  line). 
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Figure   5.7 

Change  in  absorption,  AaC^u)  on  photoexcitation  of 
polyene  with  10  sites. 
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CHAPTER  VI 
DYNAMICS  IN  DEFECTED  SYSTEMS 


It  was  noted  in  Chapter  I  that  the  electrical  conductivity  of 
trans-polyacetylene  can  be  varied  over  twelve  orders  of  magnitude  by 
adding  small  quantities  of  dopants.^-"  Clearly,  therefore,  any  model 
must  include  impurity  effects  for  a  full  understanding  of  the  transport 
properties  of  the  real  system.   Experiments  show  that  the  conductivity 
can  be  varied  reversibly  over  its  full  range  of  values  indicating  that 
the  dopant  impurities  are  non-substitutional  and  non-chemically  reactive 
with  the  0  electrons  of  the  polymer  backbone.   It  can  thus  be  assumed 
that  extrinsic  defects  only  affect  -n   electron  dynamics. 

Here  two  qualitatively  different  types  of  defect  are  modeled. 
The  site  defect,  which  may  arise  from  an  extrinsic  impurity  or  from 
morphological  defects,  such  as  cross-linking,  hybridizes  an  orbital  in 
the  defect  to  the  -t  orbital  of  a  single  carbon  atom.  The  bond  impurity 
hybridizes  an  impurity  in  the  defect  to  the  it  orbital  of  the  bond 
joining  two  neighboring  carbon  atoms.   It  can  be  thought  to  model 
extrinsic  impurities  or  cis-like  and  amorphous  regions  in  trans- 
polyacetylene. 
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Static  Model  for  the  Single  Site  Impurity 


An  impurity  of  strength  Vq  at  site  M  can  be  modeled  by  the  addition 
to  the  SSH  Hamiltonian  of  the  terra  -^ 


H.  =  V^  C^""  C,,.  (5.1) 

s    0  Ms   Ms . 


This  breaks  the  electron-hole  symmetry  of  the  SSH  equations, 
modifying  a  single  diagonal  element  of  the  hopping  matrix.  The 
locations  of  the  impurity  induced  levels  can  be  calculated  using  a 
simple  continued  fraction  scheme.   '' 

The  electronic  density  of  states  is  given  by 


p    E)  =  -  -  Lim   Im  G   5+in)  (6.2) 

^nn       7T   „      nn 


where  G^^   is  the  single  particle  Green's  funtion.   Here  the  Green's 
function  is  calculated  in  the  real  space  representation,  with  G   being 
the  propagator  for  all  paths  starting  at  and  returning  to  site  n.   From 
G^^   the  local  density  of  states  can  then  be  determined. 

As  a  simple  illustration  of  the  continued  fraction  scheme  the  band 
structure  of  the  purely  dimerized  infinite  chain  is  calculated  here. 
From  Dyson's  equation  Gqq  can  be  written  as  (see  Appendix  E) 


Gqq(z)  =  [z  -  Z^(z)  -  l^iz)]  (6.3) 


where  Z  is  the  complex  energy  E+in  and  I„(z)  and  E,  (z)  are  the  self- 

K  L 

energies  for  propagation  to  the  right  and  left  of  the  original  site 
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(0).   Now  for  t^  =  tg  +  A  /2  and  tg  =  t^  -  A  /2  then 


l^iz)    =  t^  /(z  -  Z^^(z)]  (6.4) 


where 


E^^(z)  =  tg/[z  -  ^22(2)].  (6.4) 


This  series  is  infinite.   However  as  Z„„(z)  =  Z^Cz)  on  an  infinite  line 
the  continued  fraction  has  a  periodic  form  and  is  readily  surnmable  as 


.2  r,  _  .2,,  _  .  ,,,^-^^-^ 


Zj^(z)  =  t^  [z  -  t^(z  -  Zj^(z))  ']"'.  (6.5) 


Solving  for  En(z)  gives 


'^  r-  /  \    /  2  .2    ,2,    f,  2    ,2   ,.2x2    ,,.2   2-11/2  .^  ^v 

2zEpj(z)  =  (z  +t^  -  tg)  ±  [(z  +  tg  -  tg)   -  4t^  z  J  1,0.6) 


The  appropriate  sign  of  the  discriminant  can  be  determined  from  the 
condition  that  the  density  of  states  be  positive  in  the  bands. 
Likewise 


Z^(z)  =  t^  [z  -  Ej^(z)J  '.  (6.7) 


Substituting  for  E  (z)  yields 

Li 


Gqq(z)  =  [z  -  t^  [z  -  Ej^(z))~^  -  Ej^(z)]  \  (6.8) 
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There  can  only  be  positive  contributions  to  the  density  of  states  when 
G^q(z)  is  complex.  G^q(z)  is  complex  when  E  (z)  or  5:„(z)  is  complex. 
After  a  little  algebra  this  yields  bands  in  the  regions 


E  <  2t  ,  E  >  A  ;  E  <  -A  ,  E  >  -2t  .  (5.9) 

o       o        o         o 


This  shows  the  expected  band  structure:  namely,  two  bands  separated  by 
the  band  gap  of  width  2A  .   It  is  straightforward  to  show  that  there  are 
tne  expected  square  root  singularities  in  the  density  of  states  at  the 
band  edges. 

The  single  site  impurity  is  modeled  by  adding  energy  Vq  at  a  single 
site  (taken  here   to  be  the  origin").  This  gives 


Gqq(z)  =  [z  +  Vq  -  E^(z)  -  Zj^(z)]  (6.10) 


with  Z     and  E,  as  above.  The  band  edges  have  the  same  locations  as  in 
the  pure  system,  but  the  density  of  states  at  the  band  edges  now  has  a 
square  root  convergence.   There  can  be  further  contributions  to  the 
density  of  states  at  the  poles  of  G^q(z).  After  a  little  algebra 
(Appendix  E)  one  finds  impurity  levels  at 


^±  =  TW^  ^''l  '    ^V^O^"  ^  3ign(VQ)D]      (6.11a) 


where 


D^  =  ['4tl  -    (Vq±Yq)2]2  +  l6(VQ±YQ)2t^         (6. lib) 
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and 


Y2  =  v2.tf  (6. no 


with  the  +(-)  sign  giving  the  localized  level  above  (below)  the  Fermi 
level  energy.   Thus  for  Vq  >  0  a  level  from  the  top  of  the  valence  band 
moves  up  into  the  gap,  whilst  the  level  at  the  top  of  the  conduction 
band  moves  beyond  the  band  edge  into  the  "ultra-band"  region  (i.e.  the 
ultraband  state  lies  in  the  gap  between  the  ir  band  and  the  o  band).  For 
Vq  <  0  levels  are  removed  from  the  bottom  of  the  conduction  and  valence 
bands.   It  can  be  shown  that  these  solutions  agree  with  those  of 

fin 

Baeriswyl .  ■' 

This  study  will  focus  on  donor  site  impurities  (Vq  <  0);  all 
results  are  identical  for  acceptor  sites  with  electrons  and  holes 
interchanged. 

Dynamics  In  the  Presence  of  a  Single  Site  Impurity 

The  above  analytic  calculations  do  not  treat  the  electron-phonon 
coupling  self-consistently  since  the  lattice  and  electronic  degrees  of 
freedom  are  not  allowed  to  interact  to  produce  a  fully  relaxed  lattice 
deformation  and  electronic  band  structure.  A  fully  self-consistent 
solution  can,  in  general,  only  be  found  numerically.   The  equations  of 
motion,  derived  from  the  SSH  equations,  are  integrated  and  energy  is 
removed  at  each  timestep  (by  setting  the  instantaneous  velocity  to  zero 
at  every  site  at  every  timestep)  until  the  total  energy  reaches  a 
minimum.   It  is  assumed  that  this  is  the  fully  self-consistent  solution 
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to  the  defect-modified  SSH  Hamiltonian.   At  t=0  the  lattice  is  purely 
dimerized  and  thus  this  numerical  calculation  has  same  level  of 
consistency  as  the  analytic  calculation  above.   As  expected,  there  is 
quantitative  agreement  (better  than  ]%)    between  the  analytically  and 
numerically  calculated  locations  of  the  impurity  levels.   The 
equilibrium  electronic  structure  and  lattice  deformation  are  typically 
achieved  within  0.1-0.2  psecs.  The  intragap  level  self-traps  into  the 
lattice  and  moves  closer  to  midgap  (it  still  does  not,  however,  reach 
midgap  for  a  finite  value  of  Vq).  These  localized  electronic  levels  are 
accompanied  by  a  localized  lattice  distortion,  which  consists  of  a  sharp 
step  on  the  weak  bond  neighboring  the  impurity  and  a  localized  tail  on 
the  strong  bond  neighboring  the  impurity.  With  the  sign  of  Uq  used 
here,    this  means  that  for  an  impurity  on  an  odd  site  there  is  an 
antikink-like  step,  followed  by  a  kink-like  tail  (fig.  6.1aK  Such  an 
impurity  will  be  referred  to  as  being  in  the  "antikink  topological 
sector".   If  the  impurity  is  on  an  even  site  the  impurity  is  in  the 
"kink  topological  sector".  In  this  chapter  all  site  impurities  are  in 
the  antikink  sector  unless  otherwise  stated.   Assignment  of  the 
impurity-induced  defect  as  the  bound  state  of  two  excitations  may  seem 
suprising  but  is  justified  by  the  numerical  observation  that  the 
ultraband  level  is  localized  around  the  step,  whilst  the  intragap  level 
is  localized  around  the  tail.  This  being  so,  under  certain 
circumstances,  it  should  be  possible  to  separate  the  kink  from  its  tail 
to  produce  a  topological  kink-antikink  pair  with  one  of  the  kinks 
supported,  not  by  a  midgap,  but  by  an  ultraband  level.   Indeed  this 
possibility  is  realized  and  it  has  important  consequences  for  the 
dynamics  of  site  defected  systems. 
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In  the  case  of  the  undoped  system  the  addition  of  an  extra  electron 
to  the  purely  diraerized  lattice  forms  a  polaron.   In  the  presence  of  a 
site  impurity  the  addition  of  an  electron  into  the  intragap  impurity 
level,  i.e.  a  negatively  charged  donor  impurity,  increases  the 
localization  of  the  impurity-induced  lattice  deformation.  The  electron 
then  further  self-traps  into  the  lattice  (fig.  6.1b).   However  there  is 
still  only  a  single  level  in  the  gap,  which  derives  its  parentage  from 
the  intragap  impurity  level  rather  than  from  a  polaronic  level. 

The  addition  of  two  electrons  to  the  pure  system  creates  a 
bipolaron  which  rapidly  decays  to  a  separating  kink-antikink  pair.   In 
the  presence  of  a  site  impurity  the  step  evolves  to  a  trapped  antikink 
supported  by  the  ultraband  level,  whilst  the  tail  evolves  to  a  free 
kink,  supported  by  the  intragap  impurity  level,  which  drops  to  midgap. 
This  free  kink  is  expelled  at  the  maximum  free  kink  propagation  velocity 
by  the  antikink  trapped  at  the  impurity. 

As  the  trapped  kink  is  maintained  by  a  low  lying  ultraband  level, 
optical  transitions  to  the  valence  band  are  well  above  band  edge  making 
such  a  kink  extremely  difficult  to  detect  optically.  The  bound  state 
modes  of  the  trapped  kink  may,  however,  be  observable  in  i.r. 
experiments.  As  the  trapped  kink  is  a  topological  excitation  its  bound 
state  spectrum  should  be  only  weakly  dependent  on  the  impurity  type. 
Indeed  the  i.r.  modes  usually  ascribed  to  the  pinned  kink  are  observed 
at  the  same  frequency  for  all  impurity  types.   It  should  also  be 
remembered  that  such  a  deep  lying  ultraband  level  may  hybridize  with 
the  0  bands,  the  effects  of  which  cannot  be  calculated  without  more 
detailed  band  structure  input. 
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Clearly  one  can  expect  such  a  strongly  localized  self-trapped  state 
to  qualitatively  alter  the  dynamics  of  the  system.   In  the  next  section 
we  investigate  this  further  by  looking  at  the  interaction  of  a  kink  with 
an  impurity. 

Kink  -  Site  Impurity  Interactions 

Statics 

We  first  consider  the  evolution  of  a  static  kink  and  an  impurity  on 

nearby  sites.  A  simple  continuum  calculation  has  shown  that  the  kink 

Ac 
and  impurity  evolve  independently  "^  but  we  have  seen  that  the  lattice 

defect  of  the  even  and  odd  site  impurity  are  topologically  different  so 

they  must  be  considered  separately. 

(a)  Impurity  step  and  kink  in  opposite  topological  sectors. 

The  kink  midgap  and  the  intragap  impurity  levels  evolve 
independently,  in  agreement  with  the  continuum  calculation.^ 

(b)  Impurity  step  and  kink  in  the  same  topological  sector. 

At  t=0  the  kink-like  step  is  supported  by  a  midgap  level.  The 
two  topologically  equivalent  distortions  can  be  supported  by  a 
single  localized  electronic  level.   Thus  the  intragap  impurity 
level  remerges  with  the  conduction  band,  whilst  the  kink  midgap 
level  drops  below  midgap.  The  long  time  evolution  is  sensitive  to 
initial  conditions,  particularly  the  presence  of  other  nonlinear 
excitations.   There  are  at  least  two  branches: 
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1.  The  kink  is  expelled  by  the  impurity  and  the  system  evolves 
to  an  independent  imparity  and  kink,  widely  separated  from 
each  other. 

2.  A  "trapped  kink"  forms  around  the  impurity  supported  by  an 
ultraband  state.  The  kink  intragap  state  rejoins  the 
valence  band.  The  trapped  kink  also  affects  the  continuum 
states;  each  continuum  state  now  being  localized  only  over 
the  part  of  the  chain  either  to  the  right  or  the  left  of 
the  kink. 

Dynamics 

The  above  study  of  static  kink  and  impurity  interactions  is  clearly 
not  the  whole  story;  important  new  effects  can  be  expected  in  the 
dynamics.  To  examine  these  a  single  kink  was  boosted  towards  an 
impurity  such  that  before  collision  the  kink  reached  terminal  velocity 
and  the  impurity-induced  lattice  defect  was  substantially  evolved, 
(a)   Impurity  step  and  kink  in  opposite  topological  sectors. 

At  small  |Vq|,  (-0.25  A  )  the  impurity  induced  lattice  distortion 
is  sufficiently  small  that  the  kink  can  "ride  over"  it.  At  large  |VqJ 
(-2.5  A  )  we  saw  that,  depending  on  the  velocity  the  kink  could  be 
either  transmitted  or  reflected  by  the  impurity.  We  do  not  expect  the 
outcome  of  such  a  collision  to  be  a  simple  function  of  velocity,  as  e.g. 
resonances  of  the  kink-impurity  system  and  the  presence  {or   absence)  of 
phonons  may  also  be  important.   Indeed  there  may  be  windows  of 
reflection  within  the  transmission  region,  before  a  critical  value  of 
|Vg|  above  which  the  kink  is  always  reflected.   (In  ^     dynamics   '^  it 
has  been  shown  that  for  KK  collisions,  there  are  windows  of  reflection 
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within  the  transmission  region.   Even  more  complex  dynamics  are  observed 

for  kink-breather  collisions.) 

(b)   Impurity  step  and  kink  in  same  topological  sector. 

For  small  JVq|  (-0.25  A  )  again  the  kink  can  ride  over  the 
impurity-induced  lattice  defect  (fig.  6.2a).   Here,  however,  as  the  kink 
approaches  the  impurity  the  kink  intragap  level  drops  below  midgap. 
After  the  kink  has  passed  through  it  returns  to  midgap  and  the  kink  and 
impurity  evolve  independently,  in  agreement  with  the  continuum  result. 

For  large  |Vg|  (-2.5  A  )  the  kink  initially  accelerates  to  a  high 
velocity  (greater  than  the  maximum  uniform  propagation  velocity  in  the 
free  system)  (fig.  5.2b).  The  propagation  is  smoother,  however,  with  no 
evidence  for  a  kink  tail.  At  a  large  distance  from  the  impurity  (-5a) 
the  kink  is  accelerated  to  a  very  high  velocity  (~250v  )  and  trapped  by 
the  impurity.  This  kink  appears  to  remain  trapped  for  all  times  and  is 
supported  by  the  ultraband  level  of  the  impurity.   For  a  neutral  kink 
the  midgap  level  is  singly  occupied  and  thus  does  not  fall  back  into  the 
valence  band  but,  together  with  the  intragap  impurity  level,  forms  a 
hole  polaron. 

The  trapped  kink  also  affects  the  extended  states,  each  now  being 
delocalized  only  over  the  part  of  the  chain  either  to  the  right  or  to 
the  left  of  the  impurity.   This  decreases  the  effective  chain  length, 
which  is,  as  far  as  transport  properties  are  concerned,  is  limited  by 
the  impurity  spacing  rather  than  the  physical  size  of  the  system. 

From  the  above  studies  we  conclude  that  the  ballistic  transport  of 
charge  in  kinks  through  the  system  is  unlikely  to  be  an  important 
conduction  mechanism  in  site  defected  systems.  The  above  studies  do, 
however,  contain  a  hint  of  an  alternate  mechanism.  The  trapping  of  a 
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nejtral  kink  at  an  impurity  results  in  the  production  of  a  hole  polaron 
that  might  transport  charge  by  hopping  through  the  system.  Such  hopping 
mechanisms  are  well  knovm  in  amorphous  and  disordered  systems'  and  have 
previously  been  suggested  as  the  conduction  mechanism  in  lightly  doped 
polyacetylene.-^ 


Photoexci tation  in  the  Presence  of  a  Single  Site  Impurity 

In  Chapter  IV  it  was  shown  that  breathers  are  produced  in  the 
photoexcitation  of  defect  free  trans-polyacetylene.   Clearly  if  the 
breather  mechanism  is  to  have  any  validity,  a  plausible  production 
mechanism  for  their  production  in  the  defected  system  must  be 
identified,  and  their  robustness  established.   In  this  section 
photoexcitation  in  the  presence  of  a  single  site  impurity  is  considered, 
laying  particular  stress  on  the  effects  on  the  subband  edge  optical 
absorption  spectrum.  As  representative  examples  four  simple  cases  of  a 
single  impurity  on  a  98  site  ring  are  considered.   In  each  case  the 
initial  condition  is  the  fully  self-consistent  ("relaxed")  lattice 
distortion  and  band  structure. 

(a)  Empty  intragap  impurity  level.  Photoexcitation  into  the  intragap 
level . 

Here  (fig.  6.3)  the  impurity  traps  a  negatively  charged  kink 
supported  by  the  ultraband  level.   The  intragap  impurity  level  drops  to 
midgap  forming  a  free  neutral  anti-kink  which  is  expelled  by  the  trapped 
kink.   The  singly  occupied  level  at  the  top  of  the  valence  band  and  the 
empty  level  at  the  bottom  of  the  conduction  band  move  into  the  gap  to 
form  a  hole  polaron  on  the  opposite  side  of  the  ring  from  the 
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impurity.   (This  apparent  action  at  a  distance  should  likely  be 
attributed  to  the  adiabatic  dynamics  employed  in  these  calculations,  and 
is  probably  non-physical.) 

(b)  Empty  intragap  impurity  level.   Photoexcitation  into  the  bottom  of 
the  conduction  band. 

In  this  case  a  negatively  charged  trapped  kink  and  a  positively 
charged  antikink  are  formed  in  -0.1  psecs.  The  singly  occupied  levels 
at  the  top  of  the  valence  band  and  the  top  of  the  conduction  band  move 
into  the  gap  creating  an  exciton-like  structure  far  from  the  impurity. 
The  exciton  is  dressed  by  a  breather. 

(c)  Singly  occupied  intragap  impurity  level.  Photoexcitation  into  the 
intragap  level. 

Here  (fig.  6.4)  a  negatively  charged  free  antikink  is  expelled  by 
the  negatively  charged  trapped  kink.  A  dynamic  hole  polaron  is  formed 
at  the  impurity  and  is  expelled  from  it.  There  is  no  evidence  for  an 
internal  oscillatory  mode  in  the  polaron.  This  is  an  interesting  route 
to  creating  a  dynamic  polaron  as  direct  attempts  to  boost  a  non- 
oscillatory  polaron  have  failed. 

(d)  Singly  occupied  intragap  impurity  level.   Photoexcitation  into  the 
bottom  of  the  conduction  band. 

A  negatively  charged  trapped  kink  and  a  neutral  free  antikink  form 
within  0.1  psecs.   A  dynamic  exciton  is  expelled  by  the  impurity  and  a 
breather  forms  a  long  way  from  the  impurity. 

With  these  complicated  dynamics  it  is  not  easy  to  predict  the 
effects  on  the  photoinduced  photoabsorption.  Using  the  techniques 
described  in  Chapter  III,  the  optical  absorption  of  the  system  was 
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calculated  after  it  had  evolved  for  a  substantial  period  (-0.3  psecs) 

for  the  cases  (a)  and  (o)  above. 

As  preliminary  input  the  optical  absorption  of  the  purely  dimerized 

lattice  in  the  presence  of  a  single  charged  or  neutral  site  impurity  of 

strength  -2.5  A  was  calculated.  The  presence  of  the  intragap  level 
o 

allov;s  for  subband  edge  photoabsorption.   Figures  6.5  show  the 
difference,  Aadu),  between  the  absorption  of  a  ring  with  an  impurity  and 
the  absorption  of  the  same  ring  without  an  impurity. 

For  the  neutral  impurity  there  is  a  single  subband  edge  peak, 
arising  from  transitions  from  the  valence  band  into  the  impurity  level 
(fig  6.5a").  This  is  compensated  by  an  above  band  edge  bleaching  in  the 
energy  range  2-3  A  .   For  the  charged  impurity  the  situation  is  similar 
(fig  6.5b).   Here,  however,  there  is  a  low  energy  shoulder  to  the  below 
band  edge  absorption  peak,  arising  from  transitions  from  the  impurity 
level  into  the  conduction  band.  Again  there  is  above  band  edge 
bleaching,  but  it  is  over  a  larger  energy  range.   Due  to  the  very  sharp 
lattice  distortion  around  the  step,  a  large  number  of  Fourier  modes 
(including  high  wave-vector  contributions)  are  required  for  its 
description.   In  both  cases  the  bleaching  at  the  full  band  width  is  due 
to  the  removal  of  the  square  root  singularity  in  the  density  of  states 
at  the  band  edges. 

The  optical  absorption  was  calculated  for  the  cases  (al  and  (c) 
above:  photoexcitation  into  the  impurity  level  for  both  neutral  and 
charged  impurities.  Figures  6.6  shows  the  change,  Aa((j),  in  the  optical 
absorption  on  photoexcitation.  Below  the  band  edge  the  two  cases  give 
similar  spectra:  a  strong  midgap  absorption  from  transitions  into  and 
out  of  the  free  kink  midgap  level.   The  peak  at  -A  /2  arises  from 
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transitions  into  the  polaronic  level.   The  below  band  edge  bleaching  is 
due  to  the  loss  of  transitions  into  the  intragap  impurity  level,  which 
has  now  evolved  to  midgap.   The  sharp  bleaching  above  the  band  edge 
results  from  the  loss  of  the  low  wave-vector  states  needed  to  pnodjce 
the  polaron.   For  the  neutral  impurity  there  is  further  bleaching  over  a 
wide  energy  range,  arising  from  the  loss  of  wave-vector  states  needed  to 
describe  the  trapped  kink.   (This  bleaching  is  largely  absent  in  the 
case  of  the  charged  impurity,  as  the  Fourier  descriptions  of  the  step 
and  the  trapped  kink  are  similar.) 

Photoexcitation  with  Many  Impurities 

In  each  of  the  above  experiments  the  intragap  impurity  level 
evolves  to  midgap.   In  a  general  photoexcitation  experiment  this  need 
not  be  so.   In  particular  for  a  system  with  high  impurity  density,  but 
pumped  at  a  low  optical  intensity,  insufficient  electrons  will  be 
excited  to  allow  all  the  impurity  levels  to  evolve  to  midgap.   Thus  the 
impurity  band  will  remain  largely  unaffected.  As  a  model  of  such  an 
experiment  we  photoexcited  a  single  electron  in  the  presence  of  3  and  5 
site  impurities,  placed  randomly  on  both  even  and  odd  sites. 

As  the  initial  condition  we  used  the  fully  relaxed  groundstate  in 
which,  if  the  impurities  are  widely  spaced  (>  E,   apart)  the  intragap 
levels  are  degenerate  and  if  closely  spaced  (<  5  apart)  the  impurity 
wavef unctions  overlap  to  form  an  impurity  band  of  the  finite  width.   In 
each  case  the  impurity  induced  lattice  distortions  evolve  independently, 
(a)   Empty  intragap  impurity  levels.   Photoexcitation  into  the  bottom  of 
the  conduction  band. 
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The  singly  occupied  levels  at  the  top  of  the  conduction  band  and 
the  bottom  of  the  valence  band  are  prevented  from  forming  a  KK  pair 
(with  their  associated  midgap  levels)  by  the  presence  of  the  intragap 
imparity  levels.  Instead  the  system  evolves  to  a  broadened  impurity 
band  and  an  exciton,  accompanied  by  a  high-energy  breather.   It  is 
interesting  to  note  that  all  of  the  important  dynamics  occurs  between 
the  two  most  widely  spaced  impurities,  the  rest  of  the  chain  acting  as  a 
spectator.  Here,  again,  the  effective  chain  length  is  essentially 
limited,  not  by  the  physical  size  of  the  system,  but  rather  by  the 
impurity  spacing. 

(b)   Empty  intragap  impurity  levels.   Photoexci tation  into  the  lowest 
intragap  state. 

Here  the  singly  occupied  intragap  impurity  level  drops  below 
midgap.   A  second  imparity  level  forms  a  polaronie  excitation  between 
the  two  most  widely  spaced  impurities.   The  qualitative  change  from 
impurity  levels  to  polaronie  levels  is  allowed  by  the  collapse  of  the 
kink  and  antikink  tails  of  a  pair  of  closely  spaced  impurities,  which 
form  a  pair  of  trapped  kinks.  As  in  case  (a)  there  are  a  large  number 
of  intragap  levels — two  associated  with  the  polaronie  excitation,  three 
associated  with  inert  impurities,  and  up  to  six  others  associated  with 
breathers. 

The  Single  Bond  Impurity 

In  the  remainder  of  this  chapter  an  impurity  that  couples  to  a 
single  C-C  bond  is  considered.   As  a  model  the  SSH  Hamiltonian  with  the 
addition  of  the  term 
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"b^'^^O  f^M.lS^MS  ^^MSS.IS^  (^-12) 


0  '-  M+IS   MS    MS   M+1S 

is  used,  where  Wq  is  the  impurity  strength  on  the  bond  between  sites  M 
and  M+1 .   This  alters  the  value  of  a  single  hopping  matrix  element;  e.g. 
for  the  bond  connecting  sites  0  and  1  ,  t.  -*  tp=tfl  +Wq.  Using  the 
continued  fraction  scheme  impurity  levels  are  found  at  (see  Appendix  E) 

■^.  =  5r-l'»^'B*''((tf-t|)^^.t2t^)"=]        (6., 3) 

E 

where  the  +(-)  sign  gives  the  states  above  (below)  the  Fermi  level.  The 
constant  N  takes  the  values  +1(-1)  for  the  intragap  and  ultraband 
levels.  Of  the  four  impurity  levels  predicted  by  eq.  6.13,  a  simple 
argument  shows  that  some  are  unphysical.  In  particular  at  tg=0  eq.  6.13 
predicts  that  the  location  of  the  ultraband  levels  diverges  to 
infinity.   Clearly  a  finite  impurity  strength  cannot  result  in  an 
infinite  change  in  the  energy  of  the  system.  Thus  the  ultraband  branch 
for  0  <  tp  <  tn  is  unphysical.   Consider  now  an  imparity  on  a  weak  bond 
(t.<to).  Further  weakening  this  bond  (i.e.  W  <  0)  merely  makes  the 
system  locally  more  trans-like  and  increases  the  local  gap:  therefore 
there  should  not  be  any  intragap  states.  On  the  other  hand 
strengthening  the  weak  bond  creates  a  configuration  in  which  there  are 
three  strong  bonds  in  a  row.   This  is  topologically  identical  to  a  pair 
of  kinks  on  neighboring  sites.   Thus  one  expects  intragap  states.   By 
the  same  token,  intragap  states  are  not  expected  when  a  strong  bond  is 
strengthened  but  are  expected  when  a  strong  bond  is  weakened.   Figure 
6.7  shows  the  numerically  calculated  location  of  the  intragap  levels. 


134 


This  is  in  qualitative  agreement  with  the  above  physical  arguments  and 
in  quantitative  agreement  with  eq.  6.13. 

As  in  the  case  of  the  site  impurity  the  above  calculations  do  not 
treat  the  electron-phonon  coupling  consistently.  VJe  find  that  the 
impurity  self-traps  into  the  lattice  and  at  any  value  of  ]W  |  the 
impurity  levels  become  more  localized. 

As  the  impurity  is  on  a  single  bond,  it  changes  the  local 
environment  symmetrically.  The  induced  lattice  defect  is  thus 
topologically  trivial  (fig.  6.8). 

Kink-Bond  Impurity  Interactions 

Although,  due  to  the  electron-hole  symmetry  of  the  system, 

interaction  can  only  take  place  through  the  lattice  degrees  of  freedom, 

one  can  still  expect  the  outcome  of  a  kink-impurity  collision  to  be  a 

complicated  function  of  a  number  of  parameters.  In  figs.  6.9  we  compare 

the  dynamics  of  a  single  kink  boosted  with  the  same  initial  velocity 

towards  defects  on  strong  bonds  of  strength  -  -a  /2  and  -  -  A 

o  o 

respectively.  We  see  that  for  the  weaker  impurity  the  lattice  defect  is 
sufficiently  small  the  the  kink  can  ride  over  it,  whilst  the  large 
lattice  distortion  of  the  stronger  impurity  sets  up  a  high  potential 
barrier  and  the  kink  is  reflected. 

Photoexcitation  in  the  Presence  of  a  Single  Bond  Impurity 

V/e  performed  simple  photoexcitation  experiments  with  an  impurity  of 
strength  --  -A  on  a  strong  bond.   The  initial  condition  was  the  fully 
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relaxed  system,  in  which  there  is  a  pair  of  intragap  impurity  levels. 

(a)  Nejtral  impurity.   Photoexcitation  from  the  lower  to  the  upper 
level  to  the  bottom  of  the  conduction  band. 

Here  the  electronic  occupancy  is  that  the  same  as  that  for 
photoexcitation  of  the  undefeated  system.   Although  only  -1.1  A   (less 
than  the  energy  required  to  produce  a  KK  pair)  is  added  to  the  system, 
additional  energy  is  gained  by  the  creation  of  ultraband  levels.  The 
system  is  then  able  to  create  a  KK  pair  and  a  breather  (fig.  6.1  OK 

(b)  Neutral  impurity  level:  photoexcitation  from  the  lower  impurity 
level  to  lowest  level  in  conduction  band. 

A  pair  of  oppositely  charged  polarons  is  formed  on  either  side  of 
the  impurity  (fig.  S.!!"). 

The  optical  absorption  for  the  neutral  impurity  was  calculated, 
both  before  and  after  photoexcitation  for  case  (a)  above.  Figure  6.12 
shows  the  difference,  Aa('i)'),  between  the  optical  absorption  with  and 
without  an  impurity.   There  is  a  strong  enhancement  at  -A  /4  arising 
from  transitions  between  the  impurity  levels.  There  is  also  a  weaker 
peak  at  -3A  /2  arising  from  transitions  from  the  lower  impurity  level  to 
the  conduction  band.  This  is  compensated  by  a  weak  bleaching  over  a 
wide  energy  range  above  the  band  edge. 

Figure  6.13  shows  the  difference,  Aa(a)'>,  between  the  optical 
absorption  after  and  before  photoexcitation.  There  is  a  sharp  bleaching 
at  ~A  7^4  due  to  the  loss  of  transitions  between  the  impurity  levels.   A 
strong  midgap  enhancement  is  produced  by  transitions  into  and  out  of  the 
kink  midgap  levels.   It  has  a  high  energy  shoulder  arising  from 
transitions  within  the  breather. 
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In  conclusion,  the  bond  impurity  only  results  in  intragap  levels 
for  a  limited  range  of  values  of  the  impurity  strength.   Although  the 
impurity  reduces  the  mobility  of  a  kink  the  outcome  of  a  photoexcitation 
experiment  is  largely  unaffected. 
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Figure   6.5 

Difference  between  absorption,  Aa(co),  of  98  ring  with  and  without  (a)  a 
neutral  impurity  (solid  line)  and  (b)  a  charged  impurity  (dashed  line). 
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Figure  6.6 


Change  in  absorption,  Aa(a)),  for  photoexcitatlon  of  electron 
into  impurity  level  of  neutral  impurity  (solid  line) 
charged  impurity  (dashed  linel. 


Figure  6.7 

Schematic  of  the  location  of  the  bond  impurity  induced  levels  as 
function  of  the  impurity  strength,  for  (a)  impurity  on  a  weak 
bond  and  (b)  impurity  on  a  strong  bond.  The  shaded  region  is 
the  conduction  bond  and  the  energy  zero  is  the  Fermi  level. 
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Figure   5.12 

Difference  in  absorption,  AaCw),  between  98  site  ring 
with  and  without  strong  bond  impurity  of  strength  W^  -  -A^. 
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Figure  6.13 

Change  in  absorption,  Aa((o)  on  photoexcitation  of  electron 
from  lower  to  jpper  intragap  level  strong  bond  impurity 
of  strength  W  ~  -A  . 


CHAPTER  VII 
STATICS  AND  DYNAMICS  IN  POLYYNE''^ 


A  polyyne  is  a  linear  chain  of  carbon  atoms  with  alternating  single 

72-7M 
and  triple  bonds  i.e.  iC=C^    .       It  has  been  suggested  that  the  long 

chain  polyyne,  carbyne,  is  the  stable  allotrope  of  carbon  over  the 

temperature  range  26OO-38OO  K  and  the  pressure  range  H   x  lO"'^  -  6  x  10 

atmospheres.'-'  It  has  been  further  suggested  that  polyynes  may  be  been 

seen  in  the  infrared  spectra  of  carbon  rich  stars  and  may  also  be  a 

constituent  of  interstellar  dust.'   ' 

TO    7  Q  — ft 1 

Ab  initio  Hartree-Fock  calculations'  '      have  shown  that  the 
dimerized  polyyne  structure  iC=C^_     is  energetically  favored  over  the 
cumulene  structure  4C=C=^.  .   It  is  clear  that  the  major  difference 
between  polyyne  and  polyacetylene  is  that  polyyne  has  two  TT-electron 
bands,  which  couple  instantaneously  to  the  neighboring  carbon  atoms  in 
an  identical  way.  It  is  therefore  natural  to  formulate  a  single 
electron  theory  of  polyynes  in  the  spirit  of  the  SSH  and  TLM  models  for 
polyacetylene. 


71 

Formalism' 


Within  the  framework  of  a  one  electron  tight-binding  description, 
the  TT-electron  states  of  polyyne  are  determined  by  the  linear-chain 
Hamiltonian 
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H  =  H,  -   Y.        {t   ,   C""  ,     C     +h.c.}        (7.1) 
L       ,  '  n+1 ,n  n+1 ,0  ,X  n,o,X     ' 
n  ,0  ,  A 

The  chain  consists  of  Nq  (Nq  ■*   ")  sp  -hybridized  carbon  atoms  labeled  by 

n,  (n=1 ,2, .  • .  ,Nq) .   Each  carbon  atom  has  two  degenerate,  orthogonal, 

atomic  p  orbitals  transverse  to  the  chain  axis.   The  latter  two  orbitals 

are  labeled  by  Xi\  =  \,2).      In  (7.1),  C""   ,  and  C    ,  are  fermion 

n ,  0  ,  X      n ,  a  ,  X 

operators  which  create  or  destroy,  respectively,  an  electron  with 
spin  0  in  the  Xth  orbital  of  the  nth  carbon  atom.   To  allow  for  the 
dependence  of  the  hopping  integral  t^^^^  ^   on  the  atomic  locations,  we 
adopt  the  form  used  in  polyacetylene 


t   ,    =  t„  -  a(u   -u  )  (7.2) 

n+1 ,n    0      n+1   n 


Ht  describes  the  lattice  energy  in  the  absence  of  the  overlap  of  the 
atomic  p  orbitals 


"l  =i'^  Vn4^^n(-Vr^^'-  ^7-3) 


The  linear  chain  Hamiltonian  (7.1)  describes  two  degenerate  half- 
filled  IT  bands  which  couple  to  the  instantaneous  positions  of  the  carbon 
nuclei  in  an  identical  manner.   It  follows,  therefore,  that  if  allowance 
is  made  for  the  two-fold  degeneracy  of  the  electronic  spin  and  we  treat 
the  atomic  displacements  {u  }  as  a  classical  field,  the  ground  state  of 
H  is  a  half-filled  Peierls  insulator  in  which  the  electrons  possess  an 
effective  internal  degeneracy  N  equal  to  4.   For  the  half-filled  Peierls 
insulator  ordinarily  encountered,  e.g.,  polyacetylene,  N  =  2.  Then  for 
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'/ne  uniformly  dlraerized  configuration,    u     =   ±(-1)'^Uq 


H   =    2NKu2   -      E        [t^   -    (-^'"2au^]K.^,,'n,,'   ^-^-^  ^'^ ''^ 


where  \i   is  an  index  describing  both  the  spin  index,  o,  and  the  orbital 
index,  \    and  can  take  as  the  four  values  <4   =    (t,l),  (  +  ,1),  (t,2),  {^,2] 
Now  Fourier  transforming  to  momentum  space  by 


c,  =  1-  Z  e^^"^  c^  (7.5) 

^^       /N  k        "^^ 


gives 


H  -,  =  S   [  2t-  cos  ka  c,   c,    +  4iaU-sin  ka  c,     o,   ]  •        (7.6) 
el   .  ,    0        k,M  k,M       0       k+-iT,n  k,v^ 


Hi    can  be  diagonalized  by  defining 


and  demanding 


This  yields 


wnere 


r,    =  a,  c,    +  3,  c,  (7.7a) 

k,u    k  k,y   "^k  k+Tr,u 


[^k'./el]  =  I    ^k'  ^k'.M-  ^^-^^^ 


H  ,  =  z  c,  ,   rf ,   r,  ,  (7.8) 

el   ^,      K' .u   k'.ii  k'.y 
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S.u  =  t^k  ^  Ao^in'k^]'''  =  ±\,,  (7.9a) 


with 


e^   =  2tQ  cos  Ka  (7.9b) 


and 


Aq  =  4aUQ.  (7.9c) 


The  total  ground  state  energy  is  then 


k,u 


(7.10) 


where  the  sum  is  over  occupied  levels  only.   Minimizing  this  with 
respect  to  A     gives 


1  =  ('IY^/KNq)  Z      sin^ka  E^^^  (7.11) 

k,u 


For  A/2t  <<  1 ,  which  is  the  criterion  for  the  validity  of  a  continuum 
description,  equation  7.11  then  may  be  evaluated  analytically  to  yield 

A„  =  (8t  /2.7l8)exp(-1/2X   ).  Note  that  the  dimensionless  electron- 

2 
phonon  coupling  constant  X   =  4a  /irKt^  is  twice  that  of  the  case  N  = 

ep         0 

2.   If  we  take  as  representative  values  for  polyyne  the  values  tQ=3  eV, 
K  =  68.6  eV  A   ,  and  a  =  3  eV  A~^  we  obtain  2A   =  5.0  eV.   Clearly  the 
condition  A_/2t  <<  1  is  not  well  satisfied,  implying  strong  electron- 
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phonon  interaction  for  polyyne.   However  as  we  shall  see  numerically 
2A„  =  5.88  eV,  which  indicates  that  the  weak  coupling  limit  is  not 
prohibitively  bad. 

The  Continuum  Model 

We  have  seen  in  the  previous  section  that  the  electron-phonon 
coupling  constant  is  large  (X  ~  0.4).  Although  a  continuum 
approximation  cannot  be  expected  to  give  a  quantitatively  good 
description  of  polyyne,  it  should  be  able  to  describe  qualitatively  the 
spectrum  of  nonlinear  excitations.  We  therefore  proceed  in  the  spirit 
of  the  TLM  calculation  for  polyacetylene  by  linearizing  the  electronic 
spectrum  around  the  Fermi  vector  ±kq..   '  -^  Thus 


c  =  /a  [u  (n)e      -  iv  (n)e      ]        (7.12) 


where  a  is  the  lattice  constant.  As  the  lattice  is  expected  to  dimerize 
we  expand  the  lattice  coordinate  about  ±2k^.. 


2ik„na       -2ik  na 

u^  =  IT  [A^n)e   '        +  A(n)       ]         (7.13) 
n   ^d 


where  the  factor  1/4a  is  put  in  for  future  convenience.   The  discrete 
Hamiltonian  can  now  be  linearized  to  order  a/E,.      Terms  of  higher  order 
in  a/f  are  dropped  as  are  the  rapidly  varying  terms  of  type  (-1 )'\  This 
eliminates  acoustic  phonons  from  this  model. 
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We  then  find 


H  =  -^  I  St  A^(x^  +  Z  J  dx  'F^x)[-iv„o,  I-  +  A(x)a.]v  (x)      (7.14) 
4a  \i 


where 


and 


^  (x)  =  (j  (x),  V  (x)]  (7.15) 


v^  =  2t^a  (7.16) 

e      0 


and  a.  is  the  ith  Paali  matrix.   This  is  the  TLM  equation  with  the 
1 

single  modification  that  p  can  take  on  4  values  rather  than  2  as  in 
polyacetylene.   As  we  shall  see  shortly  this  difference  is  crucial.   By 
varying  (7.14)  with  respect  to  u  (x)  and  v  (x)   we  get 


£   u^   (x)  =  -iv  |-  u^   (x)  +  A(x)v   (x)      (7.17a) 
n,\i   n,u         F  dx  n,p  n,p 


e^  V   (x)  =  +iv„  I-   v^   (x)  +  A(x)  u„   (x)     (7.17b) 
n,un,p         Fdxn,p  n.p 


where  e    is  the  nth  eigenvalue  of  the  subscript  \^ ,      Varying  (7.14) 
n ,  (J 

* 

with  respect  to  A  (x)  gives 


A(x)  =  -(4a^/K)   Z'  [u   (x)   v*   (x)  +  cc.l      (7.17c) 
n,u   "'^^     "'^ 


where  the  sum  n  is  over  occupied  levels  only.  Defining 
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f-(x)  =  u  (x)  ±  17  (x)  (7.18) 


then  one  easily  obtains 


2  3  .2-  ,      3A(x)i„+, 

'v   — o  *   ^w       -  A  (x)  ±  v^   I    'Jf-( 

F.2     k,y  rdX-'y 


[V^*'..„-''(^'±'.T7^1f*<x>  =  »-      (^-'9) 


The  groundstate  of  this  is  A(x)  =  A„(x)  with  eigenvalues  given  by 


?  2    ^  ? 

z:        =  A„  +  v-k^  (7.20) 

K,y     0     F 


where  k  is  measured  relative  to  the  zone  edge  Tr/2a  and  the  groundstate 
gap  for  the  dimerized  system  in  2A^. 


The  Gross-Neveu  Model 


Further  progress  in  determining  the  particle  spectrum  of  equations 
(7.1i4)  is  most  easily  made  by  a  consideration  of  the  Gross-Neveu  field 
theory  in  1+1  dimensions,^  '   whose  Lagrangian  is  given  by 


N 

L^'^(x)  =  I     ?''(x)  {lY  ^  -  g^,,o(x)l  ^''(x)  -  ^  a^(x)         (7.21) 

u  d  X      UN  d 

V=l  U 

Where  'i'^(x)  is  a  fermion  field,  'i^x)   =  [f^(x),  "^^{x)]   and  o(x)  is  a 

scalar  field  and  gQ«,  is  the  Gross-Neveu  coupling  constant.  The 

S'um  V  labels  the  particle  type.   Y  (jj  =  0,1)  are  the  Dirac  matrices 
(Y„  =  a^,  Y^  =  ioJ.   Varying  (7.21)  with  respect  to  H'^(x)  gives  (xq,x^ 
=  t,x) 
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^"^0  It  ^  '^1  37  "  Sgn'^^^'  ^  '{-"(x.t)  =  0.        (7.22) 


Now  looking  for  static  solutions  of  the  type 


-ie  t 
^^(x.t)  =  e   "  'F''(n;x)  (7.23) 


gives 


[e^^Q  +  i^,  33^  -  gQN  °^-^^]  *''(n;x)  =  0         (7.24) 


where  c      is  the  nth  eigenvalue.  Again  noting  that  'F^(n;x)  is  a  fermion 

function  and  that  there  are  a  pair  of  equations  in  H'^(n;x)  and  >t'^(n;x) 

for  each  value  of  v  we  see  that  if  N  =  4  (ie.  v  =  1,2,3,4)  then  (7.24) 

has  the  same  formal  structure  as  (7.17a)  and  (7.17b).   This  equivalence 

can  be  made  exact  by  making  the  transformation  "{'^(n'x)  ■*   — [  u  (n;x) 

1        2  u 

+  v^(n;x)],  4'2(n;x)  -*   -  -i[u  (n;x)  -  v  (n;x)],  go  -♦  A  and  x  ->■   VpX. 
Then  (7.22)  becomes  exactly  (7.14), 

Clearly  however,  if  the  two  models  are  equivalent  then  their 
consistency  conditions  must  match.  Here  a  complication  arises  the 
Gross-Neveu  model  is  ultraviolet  divergent  and  requires  a  cutoff 
frequency,  A,  for  renormalization.  This  may  seem  to  be  different  from 
the  TLM  model,  but  there  a  cutoff  was  introduced  in  a  natural  way  by  the 
finite  width  of  the  ir-electron  band.  The  two  consistency  conditions  can 
be  shown  to  be  equivalent. 
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Nonlinear  Excitations  in  Polyyne 

Having  shown  that  the  TLM  model  for  polyyne  and  the  statics  of  the 
semi-classical  N=  4  Gross-Neveu  model  are  equivalent  the  well  known 
results  for  the  particle  spectrum  can  be  directly  carried  over  from  the 
Gross-Neveu  model  to  polyyne.'"-^ 

There  are  kink  soliton  solutions  with  A(x)  =  A   tanh(x/^).   This 

has  a  corresDonding  eigenvalue  soectrum  with  e    =0.   As  u  can  take 

n,u 

four  different  values  the  midgap  level  can  be  occupied  by  up  to  4 
electrons.   Thus  4  x  —  =  2  states  are  removed  from  the  filled  valence 
band.   There  is  thus  a  loss  of  two  electrons  from  the  valence  band  which 
screen  two  electrons  in  the  midgap  levels.  Therefore  a  kink  can  have 
charge  +2e,  +1e,  0,  -1e,  -2e  depending  on  whether  there  are  0,1,2,3,4 
electrons  in  the  midgap  levels  respectively.   The  characteristics  of 
these  kink  states  are  shown  in  table  7.1.   It  is  important  to  note  that 
the  spin-charge  relation  of  each  of  these  kinks  is  entirely 
conventional.   This  should  be  compared  with  polyaceytlene,  in  which  the 
spin-charge  relations  are  reversed.  The  energy  of  a  kink  in  the  Gross- 
Neveu  models  is 


NA„ 

[— ^)  (7.25) 

ir 


As  topological  constraints  require  that  a  kink  (K)  and  an 

antikink  (K)  can  only  be  created  in  pairs,  the  minimum  production  energy 

of  a  KK  pair  is  2E^  =  2(4A  /tt). 
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In  analogy  with  polyacetylene  we  again  expect  there  to  be  polaron 
solutions.   From  the  Gross-Neveu  particle  spectrum  we  find 


A(x)  =  A_  -  k^v„{tanh[k^(x+x^)]  -  tanh[ k  (x-x^ ) ] }    (7.26) 
0     0F       0     9  Q     8 


where 


Kg=  r  ^  3in9  (7.27) 


^n     -1 
x^  =  —T^   tanh   [tan0/2]  (7.28) 

0   sin3        •■      ' 


and  the  angle  3  is  constrained  to  be  in  the  regime  0  S  0  <  t/2  and  is 
determined  by 

0   =   [(n+h)/4]TT/2  (7.29) 

The  energy  of  formation  of  the  polaron  state  is 


E  (0)  =  2E,  sin0.  (7.30) 

P        k 


The  polaron  eigenspectrum  has  localized  levels  at  e    =  ±e,  (k*0). 

±.M     < 

In  (7.29)  n  is  the  number  of  electrons  in  the  upper  localized  state 
which  h  is  the  number  of  holes  in  the  lower  localized  state.  The 
definition  (7.29)  above  and  the  constraint  that  0  ^  0  <  17/2  means  that  n 
and  h  can  take  values  0,1,2,3  as  long  as  (n+h)  <  4.  The  particle 
spectrum  is  easily  calculated  and  is  shown  in  taple  7.2.   (Note  the  n  = 
h  =  0  is  missing  from  this  table  as  it  simply  corresponds  to  the  ground- 
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state.)   The  principle  polaronic  excitations  are  the  polaron  (Q  = 

±e,  e  =  ±0.92A„,  5  =  1.075^,  E.  =  0.38  A^)  the  bipolaron  (Q  =  ±2e, 

c^  =  ±0.71  A  ,  C  =  l-SS^Q,  E  =  1.80  A  )  and  the  tripolaron  (Q  =  ±3e, 

e^  =  0.38  A_,  C  =  ■'•765,  S  =  0.92  A-).   Mote  that  the  polaron  also  has 

a  stable  excited  state  with  Q  =  ±e,  e   =  0.38  A.,  E,   =  1.76^,  E' 

0.92  A_  which  is  in  every  way,  but  charge,  identical  with  the 

tripolaron. 

It  is  important  to  note  that  there  is  a  neutral  polaron  (n  =  h  = 
1).  This  "polarexciton"  (Q  =  0,  £_^  =  0.92Aq,  C  =  I-ISSq,  E^  =  1.80A  ) 
can  be  considered  as  arising  from  the  self  trapping  of  both  an  electron 
and  a  hole  into  the  lattice.  Since  the  polarexciton  has  energy 


/2E  <  2  Ap^  it  is  the  lowest  energy  excitation  of  the  polymer.  Thus  one 
expects  single  electron  photoexcitation  of  polyyne  to  result  in  the 
production  of  a  polarexciton  rather  than  a  KK  pair.   The  photoexcitation 
of  two  electrons,  on  the  other  hand  gives  the  electronic  configuration 
of  a  neutral-neutral  kink-antikink  pair,  as  in  polyacetylene. 

Adiabatic  Nonlinear  Dynamics  in  Polyyne 

To  check  the  validity  of  the  above  calculations  and  particularly  to 
investigate  the  predicted  scenarios  during  photoexcitation,  the  discrete 
Hamiltonian  was  integrated  directly.  Figure  7.1  shows  the  evaluation  of 
a  98  site  system  with  196  electrons  during  the  photoexcitation  of  a 
single  electron  from  the  valence  to  the  conduction  band.  We  see  that, 
in  agreement  with  the  continuum  calculation,  a  localized  polarexciton  is 
produced.  Figure  7.2  shows  that  two  electron  photoexcitation  results  in 
the  production  of  a  KK  pair  which  rapidly  separate.   By  comparing  with 
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figure  4,1  for  the  single  electron  photoexcitation  of  polyacetylene  we 
see  that,  as  expected,  the  kink  propagation  velocity  is  higher  in 
polyyne  that  for  polyacetylene.  We  also  see  in  fig.  7.2  that  a  breather 
is  left  behind  by  the  separating  kinks. 

Using  the  formalism  described  in  Chapter  III  we  numerically 
calculated  the  optical  absorption  on  the  photoexcited  system,  in  the 
linear  response,  dipole  approximation.   As  a  preliminary  we  calculated 
the  optical  absorption  of  the  purely  dimerized  configuration.   Figure 
7.3  should  be  compared  with  fig.  3-3  which  shows  the  optical  absorption 
of  the  dimerized  configuration  of  polyactylene.   In  figure  7.3  the 
bandgap  is  2A.  =  5.83  eV  and  the  full  band  width  is  itt^  =  12  eV.  Figure 
7.ila  and  7.4b  shows  the  optical  absorption  of  the  singly  and  doubly 
photoexcited  systems  over  the  time  period  t  -  0.275  -  0.315 
psecs  (~  10-11  phonon  periods).  Figures  7.5a  and  7.5b  show  the 
difference,  AaC'jl,  between  the  absorption  after  and  before 
photoexcitation  (ie.  Fig.  7.5a  =  Fig.  7.4a  -  Fig.  7.3). 

The  optical  absorption  of  the  polarexciton  can  be  understood  by 

comparing  with  the  absorption  of  the  polaron  in  polyacetylene.  We  find 

that  for  the  polarexciton  w  /A.  -  0.63;  tnis  differs  from  the  continuum 

value  due  to  the  discreteness  of  the  system.  We  expect,  therefore  that 

the  transition  between  the  two  localized  a    levels  will  be  seen 

P 

at  -1.2  A^;  this  line  should  be  strong.   Transitions  between  the  valence 

band  and  the  lower  polaron  level  and  between  the  upper  polaron  level  and 

(2) 
the  conduction  band,  a   '  should  be  strong  and  at  -0.37  A_.   Here  we 

(1)  (2) 

find  that  a  is  stronger  than  a    ,  presumably  because  there  are  a 

larger  number  of  available  transitions  between  the  intragap  levels  of 

(?) 
the  polarexciton  than  of  the  polaron.  The  absorption  a  ~  should  be 


175 

weak  and  observable  at  -1.53  A  .   Figure  7.5a  indeed  shows  these  general 
features.   The  large  width  of  these  peaks  is  due  to  the  oscillation  of 
the  Intragap  levels  of  the  polarexciton  during  photoexcitation.   The 
above  band  edge  bleaching  is  over  a  large  energy  range  due  to  the 
narrowness  of  the  polaron — a  narrow  excitation  requires  a  large  number 
Fourier  modes,  including  those  at  high  energy  for  its  description. 

Figure  7.5b  shows  a  strong  midgap  absorption  arising  from 
transitions  into  and  out  of  the  kink  midgap  levels.   It  has  a  high 
energy  shoulder  arising  from  transitions  between  the  localized  breather 
levels.   Again  this  subband  edge  enhancement  is  compensated  by  a  strong 
above  band  edge  bleaching  over  a  wide  energy  range. 
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Table  7.1 
Characteristics  of  soliton  states  of  polyyne. 


charge        gap-state      spin         internal 

occupancy      states       degeneracy 


e  ^^  ^  S  =  1/2  4 

0  ^  ^  S  =  1  ,0  6 

-e  _^ S  =  1/2  4 

-2e  S   =   0  1 
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Table  7.2 
Characteristics  of  polaron  states  of  polyyne. 


charge     gap-state      spin      internal     stable 

occupancy      states     degeneracy   excited  state 


3e      -•-•-• S  =  1/2  4 

•  •  •  • 

2e       •  •  S  =  1  ,0  6 
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Figure   7.3 

Optical  absorption  of  a  98  site  dimerized  ring  of  polyyne  v;here 

2A  =  5.88  eV,  W  =  1 2  eV  and  the  Lorentzian  lineshape 

°  has  width  -  A  / 10. 
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Figure   7.4 

Optical  absorption,  a(u),  of  one  electron  photoexcited  polyyne 
(solid  line)  and  two  electron  photoexcited  polyyne  (dashed  line) 
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Figure  7.5 

Change  in  absorption,  Aa(w)  on  photnexcitation  of  one  electron 
(solid  line)  and  two  electrons  (dashed  line)  in  polyyne. 


CHAPTER  VIII 
CONCLUSIONS 


It  is  generally  believed  that  the  unusual  transport  and  optical 
properties  of  the  conjugated  polymers  arise  from  the  presence  of 
nonlinear  excitations  in  the  system.  ^-'•-^   Previous  work,  however,  has 
dealt  almost  exclusively  with  the  static  properties  of  the  various 
theoretical  models.   Here  we  have  undertaken  the  first  systematic  study 
of  the  dynamics  of  the  nonlinear  excitations  of  one  of  these  models. 

In  Chapter  III  we  established  the  stability  of  the  dimerized 
lattice  and  the  single  kink,  vie   further  saw  that  the  dynamics  of  a 
single  kink  can  be  well  understood  by  a  continuum  calculation.  This  and 
the  agreement  between  the  continuum  and  discrete  calculations  of  the 
optical  absorption  of  the  kink,  breather  and  polaron  show  the  relevance 
of  the  static  calculations  for  the  properties  of  the  dynamical  system. 

The  dynamical  studies  have  not  only  validated  static  calculations 
but  have  also  introduced  a  number  of  important  new  ideas.  In  particular 
the  dynamics  provide  a  unified  picture  of  the  short  time  evolution  of 
the  photoexcited  system.  The  basic  theme  appears  in  the  photoexcitation 
of  trans-polyacetylene  in  which  a  kink  and  an  antikink  are  produced; 
they  separate  at  high  speed  leaving  behind  a  coherent  optical  phonon 
package — a  breather.  Variations  on  this  theme  appear  throughout  these 
studies.  We  find  that  the  photoexcitation  of  a  polaron  also  results  in 
a  KK  pair  and  a  breather  as  does  two  electron  photoexcitation  of 
polyyne;  whilst  in  site  defected  system  one  of  the  kinks  traps  and 
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additional  nonlinear  excitations  may  be  produced.  We  also  saw  that  in 
cis-polyacetylene  strong  confinement  may  result  in  the  production  of  an 
exciton-breather  dynamical  complex. 

The  breather  has  long  been  known  to  be  an  exact  solution  to  the 
sine-Gordon  scalar  field  theory,  in  which  it  contributes  to  the 

q  ^ 

thermodynamics  of  the  system  through  its  lattice  degrees  of  freedom. 
Here  the  concept  of  the  breather  has  been  extended  to  models  with 
fermionic  degrees  of  freedom  also,  where  the  breather  contributes  to  the 
thermodynamics  of  the  system  not  only  through  its  lattice  degrees  of 
freedom  (as  in  the  calculation  of  the  dynamic  structure  factor)  but  also 
through  its  electronic  degrees  of  freedom  (as  in  the  optical  absorption 
calculations).  We  also  saw  that  the  breather  is  remarkably  long  lived, 
with  little  tendency  to  disperse  into  incoherent  phonons;  whilst  other 
studies  have  shown  that  the  breather  only  interacts  weakly  with  a 
kink.    These  suggest  that,  although  the  breather  is  only  the 
approximate  solution  of  an  effective  Lagrangian,  it  may  be  an  elementary 
excitation  of  the  discrete  system. 

At  the  outset  of  this  work  we  stressed  the  interest  aroused  in 
conjugated  polymers  by  the  remarkable  conductivity  properties  of  trans- 
polyacetylene.  Although  we  showed  in  Chapter  III  that  the  kink  is  light 
and  highly  mobile,  in  Chapter  VI  we  saw  that  disorder  may  pin  a  kink  to 
a  single  lattice  site  or  confine  it  over  a  short  length  of  chain.  This 
strongly  suggests  that  ballistic  kink  transport  is  not  the  conduction 
mechanism  in  the  defected  system.  The  extension  of  the  present  work  to 
include  more  than  one  coupled  chain  may  provide  useful  insights  into  the 
charge  conduction  mechanisms.   (Indeed,  experiments  suggest^  that  the 
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conductivity  of  the  lightly  doped  system  nay  be  via  a  variable  range 
hopping  mechanism,  similar  to  that  in  amorphous  materials.   ) 

One  of  the  other  directions  in  which  the  present  work  may  be 
extended  is  to  add  an  additional  crystal  f"ield  to  model  the  C-H  vibronic 

modes.   Although  this  may  not  be  particularly  interesting  in 

87 
polyacetylene  '  there  may  be  interesting  interactions  between  the  bond 

order  wave  and  the  charge  density  wave  in  the  third  and  quarter  filled 

Q  q 

band  systems.   The  addition  of  damping  and  A.C.  driving  to  the  equation 
of  motion  may  result  in,  not  only  spatial  and  temporal  chaos  of  the 
lattice,  but  also  chaotic  dynamics  of  the  density  of  states.   (Such  an 
observation  of  semi-classical  chaos  may  provide  useful  insights  into  the 
ill-understood  area  of  quant'ura  chaos.  ") 

As  mentioned  in  Chapters  I  and  II  this  interest  in  dynamics  is  only 
one  thrust  in  the  continuing  efforts  to  understand  the  conjugated 
polymers.   It  is  likely  that  for  a  quantitative  understanding  of  their 
thermodynamics  and  transport  properties  a  unified  theory  including 
dynamics,  electron-electron  effects,  quantum  effects  and  morphological 
effects  will  have  to  be  constructed.  Although  this  formidable  task  is 
unlikely  to  be  achieved  in  the  near  future,  Baeriswyl  and  Maki^  have 
recently  used  a  generating  functional  scheme  to  perturbatively  calculate 
the  effects  of  electron-electron  correlations  which,  it  is  hoped,  may  be 
integrated  with  the  dynamical  studies. 

One  of  the  other  paths  of  fruitful  analytic  study  may  lie  in 
collective  coordinate  schemes.  Although  a  partial  differential  equation 
has  an  infinite  number  of  degrees  of  freedom  it  may  be  possible  to 
identify  approximate  variables  e.g.  the  translational  and  boundstate 
modes  of  the  nonlinear  excitations,  with  which  soluble  dynamical 
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equations  can  be  constracted;  an  example  of  which  is  the  calculation  in 
Chapter  III  of  the  maximum  velocity  of  a  single  kink. 


APPENDIX  A 
PARAMETERS  OF  THE  SSH  AND  TLM  MODELS 


"SSH  parameters" 


"rescaled  parameters" 


t^  =  2.5  eV 


a      =  4.1  eV  A 


tQ  =  2.5  eV 


^.Q   eV  A 


-1 


K  =  21 .0  eV  ? 
M  =  13  a .  m .  u . 
a  =  1 .22  A 


K  =  17.3  eV  A 
M  =  13  a.m.u. 
a  =  1 .22  A 


u^  =  0.042  A 


A  =  0.7  eV 
o 


Uq  =  0.1  A 


A  =  1 . 9  eV 
o 


Then  jsing  the  relations 


X    =   2ct   /irt   K 


(A.I) 


to      =  /4K/M 
o 


(A. 2) 


5   =    Vo 


(A. 3) 


v^   =   2t^a 


(A. 4) 


V      =  cj   a/2 
3  o 


(A. 5) 


189 


190 


one  gets  the  parameters  of  the  TLM  Hamiltonian 


X    =   0.19 


X    =   0.3^ 


14      -1 
oJq   =    2.43    X    10        3 


u)      =    2.25    X    10 


14      -1 


7a 


C   =    2.7a 


v^   =   9.25   X    10^   ms    ^ 

4        -1 

V     =    1  .51    X   10     rns 


v„   =  9.25   X   10^  .-ns"'' 
F 

V     =    1  .37   X    lO"^   ms~^ 


APPENDIX  B 
OPTICAL  ABSORPTION  IN  THE  SSH  MODEL 


This  appendix  outlines  Horovitz's  calculation^  of  the  optical 
absorption  in  the  SSH  model. 

Formalism 

The  change  in  the  expectation  value,  i5<0(t1>,  of  an 
operator  0(t)  when  coupled  to  an  external  Hamiltonian  Hg  ,  is,  in  the 
linear  response  approximation, 

6<3(t)>  =  i  ll^ir.^<^^J    [e^"^'H^^e~^'^\  e^'"^^  3(t)  e~^^']  |*^>e^^'     (B.l) 


+ 
where  H  in  the  zeroth  order  Hamiltonian  and  n  -*■  0  . 

For  a  system  in  an  electromagnetic  field  H   =  -  J.ACtl  where  J  is 

the  current  operator  and  A(t)  is  the  vector  potential. 

Thus 


•  ft^^,^,,  I  iHt\   -iH(t'-t),  -iHt|,„  ^   nt'  ,,^,. 
<J>  =  -1   dt'<f  e    J  e        J  e    hi*  >  e    A(t')  + 

J  -00     o '  '  o 


(B.2) 


't^  ,  ^,„  I  iHt,   -iH(t-t')^   -iHt'  ,,„  ,  nt' 

dt'  <'F  e   J  e        J  e      V  >e 
-00      o '  '  o 


For  E(t)  =  E^e  ^'^^  (3.3) 
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then  A(t)  =  -A  e  '"^  (3.4) 


Inserting  a  complete  set  of  states  | 'F  >  with  eigenvalues  E^  and  then 


integrating  over  t'  one  gets 


■iLot  -itot  E 

e  e  1      o 


<J>     =    -i     E      i^-7;; r-T-    -    ^79 =; : rl    ^-     I  < 'F      I  J  I  ■■?     >|"  (3.5l 

'i(S   -E   -cj-in)        i(E   -E   -u)-iri)      loo    '      o'     '    m    '  .      -- . 


mora  mo 


Thus  the  conductivity,  a('j)  for  a  system  of  U   sites  is  given  by 


Reo(oo)    =  ^     "^^f ,    =  ^  E    |<*„|jr^>l^    [6(e   -E   -co]-  (B.6) 

E   e  m 


6(e  -E  -(jo)] 
^   o     m      ^  -' 


For  T   =   0,   E„   S   E_,  and  u)   >   0  so 


-■m   -   -^o 


Re  a{^)    =  ^  I    kf    IjI'I'   >\^  &[e  -E  -co]  (B.7) 

Nu)        '      o '     '    m    '         ^   m     o      ^ 
ra 

which  is  the  conductivity  in  the  linear  response  dipole  approximation. 

The  current  operator 

The  charge  density,  p(n),  is  given  by 


p(n)  =  eC^:^  (B.8) 


where  e  is  the  electronic  charge.  The  charge  density  is  related  to  the 
current  by  the  continuity  equation 
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3p(n)  ^ 
3t 


j(n)  -  j(n-1) 


(3.9) 


where  the  RHS  is  the  net  flux  of  current  out  of  site  n.   For  the  SSH 
Hamiltonian,  the  electronic  part  of  the  Hamiltonian,  H^i ,  is 


el 


=  -I  [t   +  a(u  -  u   J][c''C   ,  +  C""  ,C  1.       (B.IO) 
'-  o    ^  n    n+r  ■"■  n  n+1    n+1  n-" 


Now 


3C 


?  =  [c ,»: 


9t   '-  n 


(3.11) 


and  so 


3C 


9C 


e  3t     "^  at  n   ■'n  St^ 


(3.12) 


Thus 


e  at 


[t   +  dfu  -  u   J]  C"^  ,C   +  [t   +  afu   ,  -  u  ]lC*  ,C 
'■  o    ^  n    n+K-'  n+1  n   •■  o    ^  n-1    n'^ -'  n-1  n 


(B.13) 


[t   +  a  [u  -  u   J]  C*C   ,  -  [t   +  afu   ,  -  u  1  1  C'^C   . 
'■  o     '-  n    n+V-'   n  n+1    '■  o    *■  n-1    n-^ -■   n  n- 


Using 


j(n) 


m=-<» 


ap(m) 

at 


(B.14) 


then 
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j(n)  =  -ie  [t^  .  a[u^  -  .V^)][c;c^^,  -  c;^^cj       (3.15) 


which  is  eq.  3.27. 


Note  that  for 


H   =  E   V(n-m)  p(m)  p(n)  (B.16) 

ee 

nm 


[Pn'  H^J=0  (B.17) 


and  thus  (B.15)  is  correct  even  in  the  presence  of  this  electron- 
electron  interaction  term. 

The  Conductivity 

The  eigenf unction  a  has  amplitude  f  (n")  on  site  a.   Define 

a 

eigenoperators 


f  =  E  f  (n)  C  .  (B.13) 

a      a     n 
n 


Because  f  (n)  is  a  complete  basis  set,  then 
a 


C   =  E  f*(n)  f  .  (B.20) 

n      a     a 
a 


Expanding  k  >  =  n  fl\0>  (B.21) 


then 
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* 


-'^i^lV  =  "'^  ^  ^V'^^V^.l^^^'a^"^-^S^""^  '  -"^(n* ')-%(")]     (B.22) 


where  a(Q)    is  an  unoccupied  (occupied  level) 


Thus 


Re   o(oj)    =   TTB^tJ  -IT  2    Im      J^  6(£:   -e    -u))  (B.23) 

o  toN      „  '    a ,  B  3      a 


with 


M,^g   =   Z    [1    .  f-   (u^-u^^J][f^(n)f;(n.l)   -  .^^(n^DfJcn)]  (B.24) 


when  e:  and  e„  are  the  single  oarticle  energies  and 
a      8 


E^  -  E^  =  E„  -  G  (B.25) 

m    o    B    a 


Equation  8.23  and  8.24  (eq.  3.28  and  3-29)  are  again  unchanged  by  the 
addition  of  H   . 

Sum  Rule 

The  sum  rule  can  be  established  by  evaluating  a  commutator  in  two 
different  ways.  Consider 

^  X   [p(m)  [p(n).  H]]  e^^^""™\  (B.26) 

e  nm 

Now  the  commutator, 
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[A.BC]    =    [A,B]C    +    BCA.C]  (3.27a) 


and  the  ant i commutator , 


CA,3C]    =    {A,B}C   -   3{A,C}.  (B.27b) 


Then   jsing 


[C"',C    }    =   5  (B.28a) 

n     m  n.m 


{C    ,C    }    =   0  (B.23b) 

n'   m 


one  finds 


1      r-    r    /    ^      r    ^    ^   iTiT      iq(n-m") 
—  I    [pim),    [p(n),H]J    e   ^ 

e     nra 


Z-t      +au-u      J][CC      5      ,        -C      ,C      6         1 
'    '■    o       ^   n     n+V-"-    m  n     n+1  ,m         n+1   m     n.m-' 
nm 


-[t      +   a(u        -u    )][c''c      6      ,         -  C""    ,C      5         1  (B.29) 

'•    o  ^   n-1      n^-"-   m  n     n-1  ,m         n-1    m     n.m-' 


■[t      +  a[u  -u      ^^[C^C      ,6  -   C'*"C      5       ,       ] 

"^    o  n     n+1    ■"      m  n+1      n,ra  n  m     n+1  .m-" 


+  [t      +  ct[u        -u   ]][c''c      .6  -  C'^C      (5      ,      ]}   e 

'■    o  ^   n-1      n^ -"-    m  n-1      n,m         n  m     n-lim-*' 


Iq(n-m) 


4   sin^(q/2^    H   ,  .  (B.30) 

el 


Alternatively  the  commutator  can  be  evaluated  directly  from  the 
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conservation  law. 


■p(n),H]    =   i   If  =  -i[j(n)   -  j(n-1)] 


(B.31) 


which  gives 


nm 


<a|[p(m),[p(n),   h]    e^^^""'"^^  \a>   = 


nm   6  B  a 


(B.32) 


-43in^(|)   I     Z  ^4=r-  <a|j(.Ti)|3><S|j(n)la>   e''^^"  '"^    +    [q  ^  -q]. 
nm   8   ''g   "a 


And   thus 


(3.33) 


2  I 

nm 


<a|j(m)|B><8|j(n)|a> 


E    -E 
B      a 


cos[q(n-m)]   -   -<a|H   Ja>e    .  (3.3^) 


el^ 


So   for   q   =   0 


ro3 

Reo(a))dw 
J  n 


I  V    I<Q|J|8>| 


2 


B      a 


i<0|H,JO>e2 


(3.35) 


and  thus 


Reo(Lo)da)  = 
J  n 


ire 


<H   ,> 


2N  el 


(3.36) 
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which  is  eq,  3-30.  This  is  also  correct  in  the  presence  of  H^g.  So 


numerically  the  sum  rule  is  checked  by 


2t^  Z   Im  „ |^/[e  -e 
o    '  a6'   ^3  a' 


-  <H  ,  > 

el 


(3.37) 


APPENDIX  C 
THE  ANALYTIC  BREATHER 


The  Lagrangian  for  the  dimerized  half  filled  band  is 

2  '  2 
F  24A      2(jJn 

It  is  important  to  note  that  in  question  (C.I)  a  term  of  the 

•2   2 

type  A  /A   is  omitted.   This  arises  from  dynamics  generated  by  the 

valence  band  electrons.   This  is  not  included  in  the  calculated 
numerical  studies  and  will,  therefore,  be  dropped  here. 
The  equation  of  motion  from  (C.l)  is  then 


op  I  p      II 

2X  A      F         3        2 

1 2A       Wp, 

\ie   are  looking  for  a  solution  of  (C.2)  consisting  of  a  nonlinear  package 
of  optical  phonons ,  which  are  small  deviations  from  the  purely  dimerized 

lattice.  If  an  optical  phonon  has  wave-vector  ~  e  then  it  has 

2  Q? 

frequency  Wp  -  e  .  Now  use  multiple-timescale  perturbation  theory^  to 

extract  the  part  that  varies  slowly  with  both  space  and  time.  So  let 


A(x,t)  =  Aq[i  +  6(x,t)]  (C.3) 


5(x,t)    =   £:[a(X,T)    exp   ioj^t    +   A    (X,T)    exp   -iUotl    +    e   5,(x,t) 

n  n  I 
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200 
+  e^6^{x,t)    +  0(e'*)  (C.i4) 


where  the  variables  X  and  T  are 

X  =  £X  (C.5a) 

T  =  e^t.  (C.5b) 

Now  expand  (C.2)  to  order  e  .   It  is  important  to  note  in  (C.2)  that 
differentiation  is  with  respect  to  t  and  x. 

Thus  e.g.     r-  A(X,T)  exp  i-j^t  =  £A^(X,T)  exp  iw^t.         (C.6^ 

d  X  K  A  K 


It  is  Straightforward  to  show  that  substituting  eq.  C.H   into  eq.  C.2 
gives 


0   =   E°[-i  -   in   2EyAj    +    zUo) 


2A  GO 


6,    ..         .2 
+    e 


^[&,    +  -^  -  T^  5,    _    +  U^   exp    (2iaj„t)    +    IaI^ 


1  2  12      1 ,xx        2  "K    .-   -p' 


+  -  A        exp    (-2ia)^t)]    + 


r  p 

s'[6^   .  -^  -  i^  5^^^^   .  fi   {A^exp(i.^t)   -  a;   exp(-i.^t)l 

2  3 

F       ,  *  1         A  1     I     !  2 

Y2   iA^^exp(i-j^t)    +  A^^   exp(-ia3j^t )  |    -  -^  exp(3iwpjt)   -  -    |A|    Aexp   iw^^t 

*3  2 

-1|a|2a%xp   -i.j^t   -  Yexp(-3ia3j^tW    (5^    .i^fi^^^J. 


(C.7) 


f  *  111 

iA   exp   ia)„t   ■•■  A   exp  -ioj^tjl    +   0(e    1 , 
K  K 
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Now  each  order  of  the  expansion  in  e   can  be  set  to  zero  individaally 
which  gives  for  the  e   term 


Aq  =  2E^  exp  [-  j^].  (C.8) 


Now  if  we  choose 


6,(x,t)  =  7  {a  exp  21u)  t  +  A  "exp  -  2ia)  t|  -  |a|        (C.9) 


2 

and  look  at  the  e     term,  then  the  5,    term  is  of  higher  order  in  £  and 

1  ,xx  ° 

2 
can  be  formally  dropped.  Thus  the  e  terms  cancel  exactly. 

Terms  in  Op(x,t)  appears  only  linearly  to  order  e  and  thus  merely 
represents  optical  phonon  modes.   For  the  pure  breather  solution, 
therefore,  set  6p  =  0. 

Now  substituting  for  5^(x,t)  from  equation  (C.9)  and  again  noting 

that  the  term  5,    is  of  higher  order  in  e  and  can  be  dropped  then  we 
1  ,xx 

get 


2 

0  =  e^[exp(i-a)„t)[—  A_,  -  -^  A   -  ^  |a|^a]  +  c.c.].     (CIO) 
R    (jj   T    1 2   XX    3  I  1   '        J 
H 


Now  equation  (CIO)  is  the  Nonlinear  Schrddinger  Equation,  whose 
solutions  are  well  known.   In  particular  there  is  a  spatially  localized 
solution 


A(X,T)  =  b  sech  -^^  exp  -iw^b  T/3  (C.ll) 

t,  K 


then  using  relations  (C.5)  this  gives 
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A(x,t)  =  7z   sech       exp  -iu)„e^t/3  (C.12) 


where  £  =  be 


SuDStitJting  back  into  equation  (C.4)  gives 


^  /      4.  \         -^~          1-    2/2   ex  r       ,  ,,         1-2,1 

5(x,t)    =   2e   sech  — cos    [lUptd    -  -e    )j 


3' 


(C.I  3) 


1-2     ,2  2/2ex  r     -  .r.  l-2i  _-|    n/~^^ 

+  -rz      sech  — - —  [cos  2(jj_,tl  1  -  ttC  1-31  +  0(e  ). 
3  5         R  '    3 


Ilote  that  there  is  a  time  independent  part  that  reduces  the  effective 
band  gap. 

Expanding  to  appropriate  order  in  6  the  energy  density  is 

9 

2a''  '2  2    3    4    2 

E  =  -^{^  ^  ^  -  ^  -  la  -  In  «''}  (-14) 

Trv„ '   2     2    o   24    24     ' 
r  2w  j^ 

-3 
To  order  e   only  the  first  two  terms  contribute.  They  are  easily 

evaluated  by  remembering  that  energy  is  conserved  and  choosing  a 

convenient  instant  to  evaluate  the  integrals.   This  gives  the  classical 

breather  energy  as 


2/ 2  A 

0  -  r     10-2    -4  1 

E  = e   1  -  4?  ^^    O(e^)  (C.15) 

71  27 


The  Bohi — Somraerfeld  quantization  condition  is 


J  =  J   p  dA  dx  =  2mT  (C.16) 


This  elves 
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2   f  A^ 
=  — -   dx 

TTV^  J     2 


dt 


(C.17) 


V2A, 


2m 


(C.18) 


~3 
Inverting  this  to  e  and  substituting  into  (C.15)  gives 


2  2  2 

n  TT  uJp 

E  =  no)   [  1  -  -— : 

72Aq 


(C.19) 


APPENDIX  D 
CLASSICAL  DYNAMIC  STRUCTURE  FACTOR  OF  THE  BREATHER 


Define  tne  classical  dynamic  structure  factor  by 


S(q,u))  =  ^—  J  dx^dx^dt^dt^exp  [ia)(t^-t2]  -  iqlx^-x^]]       (D.l) 

(2tt) 


<exp  iqufx  t  ]  exp  -  iqu[x  t  )> 


where  a(x,t)  is  the  lattice  displacement  at  location  x  on  the  chain  and 

t  is  the  time.   The  breather  density  is  n_(a)_;T)  and  <   >  is  a  thermal 

B   B 

average.   (A  thermal  distribution  may,  however,  not  be  relevant  due  to 
the  presence  of  nonequilibrium  creation  processes.) 
Applying  the  Baker-Hausdorf  therorem 


<exp  iqu[x^t^]exp  -  iqu[x2t2)>  =  e   exp<q  u[x  t  )u[x2t2]>        (D.2) 


where  e    is  the  Debye-Waller  factor. 


Now  looking  at  one  phonon  processes  only 


exp<q  u[x^,t^]u[x2t2)>  *  q^<u[x ^ t ^ julx^t^j >         (D.3) 


But  u  (x)  =  na  +  (-1)"  u  (x)  (D.H) 


where  n  is  the  lattice  site  and  u  (x)  is  the  staggered  lattice 

n 
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displacement.  Then 


<u[x^t^]u[x2t^]>  ->  (-1)"l  (-1)"2  u^  (Xjt^)  u^  (x^t^)      (D.5) 


The  other  terms  can  be  absorbed  into  the  Debye-Waller  factor.   Now  in 

,,  ..  ■,  ■    ■  J.    /    ,  \n  iiTx/a 

the  continuum  limit    (-1)      =  e  ,30 


(2ti)' 


S(q,a))    =   e —  J  dx   dx   dt   dt   exp[i,ij(t^-t^)-i(q-Tr/aMx   -X    )]q 


<u(x^t^ )j(x2t2)> 


(D.6) 


Assuming  that  the  breather  behaves  as  a  relativistic  particle,  i.e., 


x  -*   Y(x-vt)  (D.7) 

t  ^  yft-vx/c"] 

where  T  is  the  Lorentz  factor  and  c^  is  the  velocity  scale.  We  also 
assume  that  the  thermal  average  can  be  implemented  by  integrating  over  a 
classical  velocity  distribution,  P(v1. 


Thus 


-?W   '^R^^R''^'^      2r  ?  ?  2 

S(q,a))    =   e —  q  JdvP(v)[Jdx   dt   exp[-iY    ( q+Tr/a-vcj/c  ")x   exp   iY 


(u)-qv)t^]u(x^t^)][x^    ^   x^,    t^    >   t^]  (D.8) 


For  simplicity  assume  the  Maxwell-Boltzmann  distribution  can  be  used  so. 
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T  ->■    1  (D.9) 


P(v)  .  P^'^v)  3  [^)^/2  exp[-av2/2c2]  (D.10) 

o 


where  a  =  E  /kT. 


Now  for  the  breather 


u(x,t)  =  2c(a)_)  sech(x/d)  cos  w^t  +  ^  e  (u)„)  sech  (x/d)[cos  2u)„t  -  s] 

d  □       J        D  D 


wnere 


(D.n) 


1  ~2^ 


and  (D.12) 


d  =  ^/2/2  I. 


Define  the  q-dependent  from  factors 


f\q)    =  Jdye   ^^^3ech(y/d)    =   ird  sech(TTqd/2)  (D.13) 


f^(q^    =   jdye   ^^^seGh^(y/d)    =   irqd^co3ech(TTqd/2)  (D.14) 


207 
Now  performing  the  integrates  over  x,  and  t-.    in    (D.8)  gives 


S(q,a3)  =  ^  ^'    e  ^'"'jdvP™(v)  {if^l<)[5(a3+a)  -qv)  + 

6(a)-a)g-qv)]  +  -  l^r    (k)  [-6((u+2ojg-qv)  +  -  &{oi-2ui^qv)  (D.15) 


-35  (a)-qv)]J  dx^dt  exp[i(q+Tr/a-va)/G^^x  ]exp[i(u)-qv)t  ]u(x  t  ) 


with  k  =  q+Tr/a  -  vw/c''  (D.lb) 


Now  perform  the  v  integral  followed  by  the  integrals  over  X2  and  t-,  (the 
integral  over  tp  introduces  a  5-function  which  forbids  terms  coupling 
e.g.  uj+oj  with  (j-u)  ) . 
One  finds 


S(q,.)  =  e  ^'     ^     \   [i'{|f^(k]|^  ?'"(--^)  - 
(2Tr)^  '  '^ 


(D.17) 

|f'[k-]|^  P''^"[^]  -  36  |f^[kJ|V"(./q)l] 


,2r  -1  |2  ^NRr'^  ^"b-,    ,^  i^2r,  1  1  2  I>JR , 


-1  ij±nu 

where  k*  =  d   [q  +  ir/a  -  ^  [ -]]  (D.181 

n       '■  ^  2  ^   q   ^  -^ 

o 

_l  2 

k  =  d   [q  +  u/a  -  -^).  (D.19) 

o  2 

c  q 
o  ^ 


APPENDIX  E 
CONTINUED  FRACTION  SCHEME  FOR  IMPURITY  LEVELS 


Formalism 


The  local  electronic  density  of  states  at  site  n  is 


p  (E)  =  —  Lim  Im  G   (E+in)  (E.I) 

n-^O 


where  Gj^^^  is  the  single  particle  Green's  function.   For  the  Harailtonian 


H=Ze   C   C   +St     ,(C   ,C   +h.G.)  (E.2) 

n  n  n      n,n+1   n+1   n 
n  n 


Dyson's  equation  may  be  written 


G(z)  =  (z-H)  (E.3) 


=  1  +  ^  HG(z)  (E.n) 

z   z 


Taking  matrix  elements  between  the  basis  of  Wannier  functions  {|i>} 
gives 


z  <  i|G(z)|j>  =  <i|j>  +  G^<i|G(z)|j>  +  E  <i|t^  ^i+l  '  ^'^'^'^ '^'-^^       ^^■^• 

k      ' 
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which  may  be  written 


zG..(z)  =5..  *  ..  a..(z)  .  E  t.^,  G,^.(z)  (E.6^ 


Iterative  elimination  of  off-diagonal  terms  gives 


where 


°-'^' =  ^-^-^T^ 


t2 


L  ?  (E.8a) 

z  -  e .  ,    t .  ,  .  „ 
1-1  -   1-1 ,1-2 


and 


^2 


R     _      _  .2  (E.8b) 

^   "-i  +  l   ^i4-l,i  +  2 


Undefected  Lattice 

Let  us  calculate  the  density  of  states  at  site  zero  for  the 
undefected  lattice.   Let  t2|  2i  +  l  "^  ^A '  ''21  2i-1  "  ''B 

then 


Z  (z)  =  ^ (E.9) 

^      z-Z^^(z) 
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and 

E  "(z)  =  —^ .  (E.IO) 

Since  E„„(z)  =  E„(z)  on  an  infinite  line  then  the  continued  fraction  is 
periodic  and  readily  summable 


E^(z)  =  t2  [z  -  t^  [z  -  l^(z)y^r\  (E.11) 


Solving  for  E„(z)  gives 


2zE^(z)  =  (z^  +  tj  -  t^)  ±  [(z^  +  t^  -  t^)2  -  nt^z^]^^^     (E.12) 


Using 


gives 


Gqq{z)    =  [z  -  t2[z  -  Sj^(z)]~^  -  Ej^(z)]    .         (E.U) 


From  equation  (E.I)  there  are  contributions  to  the  density  of  states  if 
Gqq(z)  is  complex.  Combining  (E.12)  and  (E.14)  yields  contributions  to 
the  density  of  states  at 


(z^  +  tl   -   tf)^  -  '4t^z^  <  0  (E.15a) 

D       A  D 
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2     ?     7    7  5  ? 

{tT   -  tg  +  tp^  -  4t-z'^  <  0  (E.15b) 


Assuming  t»  >  tg  this  gives 


E  <  t^  +  tg  (E.l6a) 


E  >  t^  -  tg  (E.l6b) 


S  <  -t.  +  t„  (E.16g1 

A       D 


E  >  -t^  -  tg  (S.iDd) 


Now  for  t„  =  t„  +  A  /2  and  to  =  t_  -  A  /2  the  bands  are  in  the  region 

A      O      o  d      O      o 


A   <  E  <  2t    and   -A   >  E  >  -2t  .  (E.l?) 

o        o  o         o 


That  is  we  have  a  full  band  width  of  W  =  ^t^  and  a  band  gap  of  2A 


The  Site  Defected  System 


For  an  impurity  at  site  0  we  add  to  the  SSH  Hamiltonian  the  term 


H  =  -V-  c!;  C^  (E.18) 

s     0  Os  Os 


then 


%o^^)  =  rTi^-^Tjir::-ViT  ^^-^9) 
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where  E,  (z)  and  E_(z)  are  given  by  equations  E.I  2  and  E.13  above.   There 


:an  be  contribJtions  to  p^(E)  at  the  poles  of  Gqq(E) 


o 


z  +  Vq  -  Zj^(z)  -  Ej^Cz)  =  0  (E.20) 


It  is  easy  to  show  that  this  jives 


wnere 


d2  =  [^^2  _  ^^^  ^  ^^^2^2  ^  ^g^^^  ^  y^)2^2        ^g^21b) 


and  Yq  =  'it^  +  Vq  (E.21C) 


where  L  and  M  can  take  the  values  of  ±1 . 
Expanding  eq.  E.21  to  lowest  order  gives 


E  =  L(t^  +  Mtg)(l  -  V^/t^)  (E.22) 


This  identifies  the  various  band  edges  v/ith  particular  signs  of  L  and  M, 
and  shows  that  for  Y^  <  0  states  are  removed  from  the  bottom  edges  of 
the  conduction  and  valence  bands,  whilst  for  V^  >  0  state  are  removed 
from  the  top  edges  of  the  valence  and  conduction  bands. 
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Bond  Impurity 

The  bond  impurity  modifies  the  transfer  integral  between  two  single 
sites  e.g.  site  0  and  site  1. 


So 


t„,  =  t^  =  t„  +  W  (E.23) 

01     E    A    o 


A  little  algebra  then  yields 


E 


Using  the  physical  arguments  of  Chapter  VI  we  find 


[tgl  <  t^  ,  t,^  >  tg  (E.25a) 


Itgl  >  t^  .  tft  <  tg  (E.25b) 


and 


N  =  +1  |tg|  >  t^  (S.26) 
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